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ABSTRACT 

We present a database of structural and dynamical properties for 153 spatially resolved star clusters in the 
Milky Way, the Large and Small Magellanic Clouds, and th e Fornax dwa rf s pheroidal. This database c omple- 
ments and extends others in the literature, such as those of Harris ( 1996) and Mackev & Gilmore (200 3albich . 
Our cluster sample comprises 50 "young massive clusters" in the LMC and SMC, and 103 old globular clusters 
between the four galaxies. The parameters we list include central and half-light averaged surface brightnesses 
and mass densities; core and effective radii; central potentials, concentration parameters, and tidal radii; pre- 
dicted central velocity dispersions and escape velocities; total luminosities, masses, and binding energies; 
central phase-space densities; half-mass relaxation times; and "«-space" parameters. We use publicly available 
population-synthesis models to compute stellar-population properties (intrinsic B — V colors, reddenings, and 
V-band mass-to-light ratios) for the same 153 clusters plus another 63 globulars in the Milky Way. We also take 
velocity-dispersion measurements from the literature for a subset of 57 (mostly old) clusters to derive dynami- 
cal mass-to-light ratios for them, showing that these compare very well to the population-synthesis predictions. 
The combined dataset is intended to serve as the basis for future investigations of structural correlations and the 
fundamental plane of massive star clusters, including especially comparisons between the systemic properties 
of young and old clusters. 

The structural and dynamical parameters are derived from fitting three different models — the modified 
isothermal sphere of|King| l|1966): an alternate modified isothermal sphere based on the ad hoc stellar dis- 
tribution function of Wilson (1975); and asymptotic power-law models with constant-density cores — to the 
surface-brightness profile of each cluster. Surface-brightnes s data for the LMC, SMC, and Fornax clusters are 
based in large part on the work of Mackev & Gilmore (2OO3g03), but include significant supplementary data 
culled from the literature and important corrections to Mac key & Gilmore's V-band magnitude scale. The pro- 
files of Galactic globular clusters are taken from Trager, King, & Diorgovski ( 1995). We address the question 
of which model fits each cluster best, finding in the majority of cases that the Wilson models — which are spa- 
tially more extended than King models but still include a finite, "tidal" cut-off in density — fit clusters of any 
age, in any galaxy, as well as or better than King models. Untruncated, asymptotic power laws often fit about 
as well as Wilson models but can be significantly worse. We argue that the extended halos known to character- 
ize many Magellanic Cloud clusters may be examples of the generic envelope structure of self-gravitating star 
clusters, not just transient features associated strictly with young age. 
Subject headings: globular clusters: general — galaxies: star clusters — Magellanic Clouds 



1. INTRODUCTION 

Next to elliptical galaxies, globular clusters are the 
most thoroughly modeled and best understood class of 
"hot" stellar system. The simple models of single- mass, 
isotrop ic, lowered isothermal spheres developed by iKing l 
(1966) have been fit to a large majority of the ~ 150 
Galactic globulars currently known, yielding comprehen- 
sive catalogues of clust er structural p aramet ers and de- 
rived physical propertie s llDiorgovskilll993tlPryor & Mevlan 
ll993lrTrager. King. & Pi orgovskil995tlHarrisl 19961) . These 
have been used to explore a multitude of scaling rela- 
tions and in terdependences between the various properties 
( Diorg ovski & Mevlanl 11994). leading to the definition of 
a fundamental plane for globular clusters which is anal- 
ogous to but physically distinct from that for early-type 
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galax i es and bulges (e . g.. lDiorgovskilll995t iBurstein et alJ 
Il997t Bellazzinil 1l998t iMcLaughlinl 120001) . Spectroscopic 
and HST imaging data have also been collected, and fit 
wilh lKingl dl966l) models, fo r a la rge number of globul ars in 
M31 (e.g..lDiorgovsMetai1[1997tlDubath & GriUmairll 19971 
iBarmbv. Holland & Huchral2002l) , M33 llLarsen et alJ2002l). 
and NGC 5128 = Ce ntaurus A dHolland Cote. & Hessed 
U999tlHarris et alJl2002HMartini & Holl2004l) . These clusters 
appear to follow essentially the same scaling relations and lie 
on the same fundamental plane as Galactic globulars. 

That these data contain important information on the for- 
mation and evolution of globular clusters (GCs) is clear. But 
cleanly separating formative from evolutionary influences on 
the present-day form of GC structural correlations and the 
fundamental plane is difficult in the absence of any rigor- 
ous, first-principles predictions for the systemic properties of 
newly born star clusters. This issue only gains in importance 
with the emergent consensus that the "super" star clusters now 
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forming in nearby starbursts and galaxy mergers, and even 
massive young clusters in more quiescent disks, may well be 
close analogues to what globular clusters were a Hubble time 
ago. What, then, can the properties of GCs tell us about how 
star clusters form in general? Conversely, what can a compar- 
ison between the properties of young clusters and GCs add to 
our understanding of the globulars? 

Important steps towards addressing these questions have 
been made by studies of the mass-radius re lations for young 
massive clusters in merging galaxi es (e.g.. |Zepf et alJl!999l) 
and relatively nearby spiral disks (Larsen 2004), and, on a 
somewhat different s cale, for the (young) nuclear clusters in 
very late-type spirals (Bo ker et alJ2004lWalcher et alJ' 2005). 
However, the data on these clusters — all of which are at dis- 
tances of many Mpc — are necessarily much less complete, 
and in some respects more uncertain, than those for globu- 
lar clusters in the Milky Way and even out to NGC 5128. 
A comprehensive study of young cluster properties vs. GC 
properties requires a catalogue for the former which includes 
the full suite of physical parameters routinely calculated for 
Galactic GCs, preferably obtained within the same modeling 
framework and to comparable precision and accuracy. 

Our main goal in this paper is to contribute to such a 
database, focusing specifically on the nearby populations of 
massive clusters in the Large and Small Magellanic Clouds. 
In particular, we fit a variety of models to obtain derived 
parameters for 53 LMC clusters and 10 SMC clusters with 
high-quality, HST-based surf ace-density profiles of t heir in- 
ner parts measured by Mackev & Gilmorel ll2003albl) . Thir- 
teen of these objects are globular clusters older than 10 Gyr; 
the other 50 ran ge in age from several Myr to a few Gyr. 
At the same time, Mackev & Gilmore (2003c) have published 
new surface-brightness profiles for the 5 old GCs in the For- 
nax dwarf spheroidal, so we include these in our catalogue 
as well. As we describe below, some aspects of our analyses 
of these LMC, SMC, and Fornax clusters go beyond what is 
normally done for GCs. Thus, we also take surface-brightness 
profiles for globulars in the Milky Way from the collection of 
Trage rTKing. & Diorgovskil (119951). and some supplementary 
data from the catalogue of lHarrisH l996). in order to model as 
many Galactic GCs as is practical in precisely the same way 
that we treat the Mackey & Gilmore datasets. The exact num- 
ber of Galactic GCs that we include depends on the context of 
our various calculations. 

Much in the existing catalogues of Galactic GC prop- 
erties is based on the same sur face-brightness data from 
Trager. King. & Diorgovski ( 1995) that we employ here; and, 
as we said above, the catalogues already contain the results 
of reasonably uniform modeling of the clusters. However, all 
of this is within the context of the theoretical King (1966) 
isotropic, m odified isothermal spheres. Mackev & Gilmore 
(2003a b c), for their part, tabulate estimates of cluster cen- 
tral surface brightnesses, core radii, total luminosities, and 
total masses derived from fits of power-law models for their 
surface-density profiles: I(R) oc [1 + (7?/V ) 2 r (7 ~ 1) < /2 . (Mackey 
& Gilmore's cluster masses further depend on their applica- 
tion of population-synthesis models.) One of our aims here 
is to homogenize the analysis of the basic data on young and 
globular clusters alike, and thus we completely re-fit all of the 
LMC, SMC, Fornax, and Milky Way cluster profiles with both 
iKing l ill 9661) and power-law models. It is then important to 
recognize Mackey & Gilmore's motivation for fitting asymp- 
totic power laws to their cluster sample in the first place. 



In their pionee ring study of the structure of young clus- 
ters in the LMC, lElson. Fall & Freemar] dl987l) found that 
[KiTT^ dT966l) models are not always capable of describing the 
outer-envelope structures of these objects, which are very ex- 
tended spatially and may be better fit by power laws that do 
not include the sharp tidal cut-offs built into King models. 
Some subsequent studies of other yo ung massive c l usters have 
come to the same conclusion (e.g., Larse r] l2TMISchweIzeil 
2004). It has been argued that this reflects the presence of 
rather massive halos of unbound stars around the clusters, 
which (it is presumed) will eventually be stripped away by 
tides to leave behind a more familiar, King-like body. How- 
ever, so far as we are aware, there have not been any sys- 
tematic attempts made to fit large samples of young clus- 
ters w ith structural models that are more extended than lKingl 
( 1966) but still spatially truncated rather than formally infi- 
nite like power laws. Nor, for that matter, has any such test 
been performed on old GCs. Thus, in this paper we also fit 
our combined LMC/SMC/Fornax/Milky Way cluster sample 
with a third type of model which is intermediate to King and 
power-law models. For this we have chosen a spheri cal and 
isotrop ic version of the model originally developed bv lWilsonl 
( 1975) for application to elliptical galaxies. Wilson's model 
is essentially a single-mass, modified isothermal sphere like a 
iKind i ll 9661) model in its core but with an ad hoc change in 
the "lowering" term in the stellar distribution function to give 
more extended (though still finite) halos for otherwise similar 
clusters. 

In ^2 below, we discuss in detail the surface-brightness 
data that we use to fit models to the clusters in our sample. 
In the LMC and SMC, we supplement many o f the HST- 
derived profiles in Mackev & Gilmore (2003a b) with older 
(ground-based) starcount densities for the same clusters. In 
many cases this allows the density profiles to be defined out 
to substantially larger proje cted radii. For the Fornax cl usters, 
we rely exclusively on the Mackev & Gilmore ( 2003cp data. 
We work in the V band, and in all cases we check Mackey 
& Gilmore's magnitude scales against ground-based aperture 
magnitudes from the literature. We find it necessary to re- 
calibrate Mackey & Gilmore's surface brightnesses, since as 
published they lead to integrated magnitudes that are always 
too faint — often by several tenths of a mag — relative to the 
ground-based numbers. The zeropoint correction for every 
cluster is tabulated. Our discussion of the Milky Way GC 
surface-brightness data is very brief, as we simply take them 
(also in the V band) from Tra ger. King. & Diorgovskil ( 119951) 
with minor modifications. 

Following that, in ^3] we explore the use of population- 
synthesis models to infer an extinction Ay and mass-to-light 
ratio Ty for each of the 6 8 LMC, S MC, and Fornax clusters 
from Mackev & Gilmore (2003a be) — who have also pro- 
vided an age and metallicity for every cluster, either compiled 
from CMD studies in the literature or estimated themselves — 
and for 148 Galac tic GC s with metallicities listed in the cat- 
alogue of Harris ( 1996). We resort to the use of extinc- 
tions based on population-synthesis model colors to derive 
intrinsic luminosities, etc., from the observations of LMC 
and SMC clusters, because direct measurements of redden- 
ing are not available for the majority of that sample. (Mea- 
sured extinctions are available for the Fornax and Milky 
Way GCs, and we use those as a check on our population- 
synthesis model values.) We need theoretical mass-to-light 
ratios for most of the clusters in all four galaxies in or- 
der to convert from luminosity to mass when determining 
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a number of dynamical cluster parameters (velocity disper- 
sions have been observed for only a fraction of the clusters 
under consideration here). We tabulate Ay and Yy for ev- 
ery cluster as determined using each of two publicly avail- 
able population- synthesis code s llFioc & Rocca-Volmerang el 
1997; iBruzual & Charlotll2003l) under a number of different 
assumptions on the form of the stellar IMF. For our subse- 
quent modelin g, however, we only use th e values predicted 
by the code of IBruzual & Charlotl (2003) with the disk-star 
IMF of lChabrierl i )2003l) . iMackev & Gilmorel J2003albld) also 
used population-synthesis modeling to estimate mass-to-light 
ratios for their LMC/SMC/Fornax clusters; but they assumed 
a much steeper IMF than we do and obtained systematically 
different numbers for Ty. 

Section[4]first describes (" 34.lt some sal ient aspects of King, 
power-law, and Wilson models. Then 34.21 shows in detail 
how each compares to the observed surface-brightness and 
internal velocity-d ispersion profiles of the well s tudied GC ui 
Centauri (see also McLaughlin & Mevlan 2003), before pre- 
senting the bulk of our results, the basic parameters of each 
model fit to each of the 68 LMC, SMC, and Fornax clus- 
ters and 85 Galactic globulars. These are accompanied by 
a number of derived cluster properties, all evaluated within 
each of the three fitted models: central and half-light aver- 
aged surface brightnesses and mass densities; core and effec- 
tive radii; concentrations and tidal radii (or infinite power-law 
slope); predicted central velocity dispersions and escape ve- 
locities; total luminosities, masses, and binding energies; cen- 
t ral p hase-space densities; and half-mass relaxation times. In 
34.3l we compare our fits to those in existing catalogues (when 
the latter overlap with our work), and we check that our de- 
rived cluster properties are generally well defined regardless 
of which model is fit to the data . 

We also investigate in 34.31 the question of which model 
tends to fit these clusters best. We find that, for w 90% of our 
full sample of young massive cluster s and old globular clus- 
ters, the more extended Wilson ( 1975) models provide equally 
good or significantly better fits than King ( 1966) models. Un- 
truncated power laws generally describe the young clusters 
about as well as the spatially limited Wilson models; only oc- 
casionally are they slightly better. Thus, we conclude that (1) 
unlimited power laws are not the only description possible for 
the halos of young LMC/SMC clusters in particu lar, and (2) 
structure extending beyond what is predicted by King ( 1966) 
models is a fairly generic feature of any star cluster, rather 
than something present only at young ages. It is not obvious 
that all such extended halos must necessarily be unbound, but 
it is a subtle problem to measure accurately the tidal radius of 
a cluster and to compare it properly to the theoretical Roche 
lobe defined by a galaxy's tidal field. We do not tackle the 
issue in any detail in this paper. 

In |3 we collect velocity-dispersion data from the litera- 
ture for a subset of 19 Magellanic Cloud and Fornax clusters 
and 38 Galactic GCs, and use these in conjunction with our 
surface-brightness fits to derive dynamical mass-to-light ra- 
tios for them. We compare these directly to the population- 
synthesis model values which we used to derive all mass- 
dependent cluster properties. The agreement is very good on 
average. 

Finally, 50 presents the "«-sp ace" paramet ers of 

iBender. Burstein. & Fabeil i 19921) and iBurstein et all i 19971) 
(or, more precisely, mass-equivalent versions of these) for 
the 153 young and globular clusters to which we fit all three 
of our structural models. We also tabulate in this Section 



the galactocentric radii of all the objects. These quantities 
complete what is needed to construct the fundamental plane 
of star clusters in any of the equivalent formulations that can 
be found in the current literature. However, we leave the ac- 
tual delineation and interpretation of all cluster correlations, 
including discussion of young vs. old populations, for future 
work. 

2. STARCOUNT DATA, APERTURE PHOTOMETRY, AND 
SURFACE-BRIGHTNESS PROFILES 

In this Section we collect and combine, in a uniform fash- 
ion, available data defining the run of surface brightness as a 
function of radius in a significant number of nearby and well- 
resolved globular and young, massive star clusters. 

In the case of the Milky Way globular cluster (GC) 
system, this task has in fact alr eady been completed by 
Trager, King, & Diorgovski ( 1995), who gathered inhomoge- 
neous surface-brightness (SB) and starcount data from the lit- 
erature for 124 GCs and combined them objectively into a 
single, zeropointed V-band SB profile for each cluster. 1 This 
database has served as the raw materia l for the standard cat- 
alogues (D iorgovskill993HHarrisl 19961 updated 2003 Febru- 
ary 2 ) of King ( 1966) model par ameters for Galactic GCs. 
Here w e similarly work from the Trager . King. & Diorg ovski 
i 19951) GC profile data for our modeling of Milky Way glob- 
ulars. Although some of these profiles have been superseded 
by more recent work on a few individual clusters, we have 
not attempted to incorporate any such updates into the Trager 
et al. catalogue, preferring instead to draw on data that have 
been processed in a homogeneous way by a single set of au- 
thors. A few further deta ils on our handling of the Trager et 
al. data are given in 32.2l below. 

The situation for well-resolved massive clusters in the 
Magellanic Clouds and the Fornax dwarf spheroidal is not 
quite as simple, although in th ese cases the recent studies of 
Mackev & Gilmore (2003 333) h ave gone a long way towards 
provid ing a close analogue to the Trage r. King. & Diorgovski 
( 1995) dataset for Milky Way globulars. Mackey & Gilmore 
performed starcounts on the inner w 100" of 68 massive clus- 
ters (including 18 old globulars) from archival HST (WFPC2) 
images, providing the most accurate definitions of the core 
structures of these objects. (In the LMC, at a distance of 
D = 50.1 kpc, 100" corresponds to 24.2 pc; in the SMC, with 
D = 60.0 kpc, 100" = 29.1 pc; for the Fornax dwarf, D = 137 
kpc and 100" = 66.4 pc.) Mackey & Gilmore further con- 
verted their starcount densities to V-band fluxes and published 
their results as standard surface-brightness profiles. 

Our modeling of the LMC/SMC/Forna x clusters is there- 
fore based in largest part on the work of Ma ckev & G ilmore 
(2003a b c). Again in the interest of confining our attention to 
as uniform a raw dataset as possible, we have not brought into 
the sample any other Magellanic Cloud clusters studied by 
other authors. (Note however, that Mackey & Gilmore's clus- 
ters do include many previously observed from the ground 
and modeled by other authors; see their papers for details of 
the overlap and comparisons with earlier work.) 

We did, however, find it useful to go back through the lit- 
erature to supplement Mackey & Gilmore's HST data with 

1 We note in passing that Trager et al . give the nu mber of clusters in their 
catalogue as 125. However, according to Harris 1 1996), the GCs Terzan 5 and 
Terzan 11, which appear separately in Traaer. Kine. & Dioreovski 1 1995), 
are in fact the same object. 

2 Available online at 

|http : / /physwww .mcmaster . ca/| fe7Eharris/mwgc . dat . 
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(1) any available ground-based starcount data at projected 
radii R > 100", and (2) all available concentric-aperture pho- 
tometry for these 68 clusters. The detailed comparison of 
the HST data against ground-based aperture photometry and 
large-radius starcounts, and our combination of these into a 
single, properly zeropointed V-band SB profile for each clus- 
ter, are what we describe now. 

2.1. LMC, SMC, and Fornax Clusters 

Table ^ lists a U sources of star counts that we have 
taken from the literature for the 53 massive clusters (in- 
cluding 12 old globulars) studied by Mackev & Gilmore 
(2003a) in the Large Magel lanic Cloud; the 1 (including 
one GC) studied by [Mackey & Gilmord (l2003bl) in the SMC; 
and the 5 globular clusters of the Fornax dwarf spheroidal 
(Mack ev & Gilmorell2003cl) . We note again that, although 
Mackey & Gilmore have published their data as SB pro- 
files, they are in fact number densities converted to net mag 
arcsec -2 ; the other papers listed in Table [2 publish their re- 
sults directly as number densities N per unit area on the sky, 
either already corrected for or accompanied by an estimate 
of background contamination. The older data can be com- 
bined with those of Mackey & Gilmore by transforming the 
ground-based N as n = C -2.5 log N, with the constant C 
chosen to provide the best average agreement with the HST 
surface brightnesses at clustercentric radii where the observa- 
tions overlap. 

Fo r each of their 68 clusters in th ese three Galactic satel- 
lites, iMacke^^^iimo^ <|2Q03ilIE hereafter referred to col- 
lectively as MG03) have produced a "primary" surface bright- 
ness profile in the V band and a "secondary" profile in one of 
B or /. We have chosen to work on the V magnitude scale, 
but to make maximal use of the new HST data, we have com- 
bined the data in the secondary bands with the main V-band 
data by applying a constant color term to the former. These 
instrumental colors are determined simply from integrating to 
find the total magnitude of every MG03 cluster in both V and 
B or /, and then adding the difference ( V - B) or ( V - 1) to ev- 
ery point in the secondary SB profile. This produces a single 
profile with two independent estimates of fiy — both used in 
our model fitting below — at every clustercentric radius. The 
instrumental colors we applied are listed for reference in Ta- 
ble [2] and are in agreement with the total cluster magnitudes 
already calculated by MG03. 

After combining the primary V and secondary B/I SB 
data for each cluster, we parsed them to be left only with 
strictly independent surface brightnesses. MG03 performed 
starcounts in four different series of concentric apertures cov- 
ering each cluster, each series using a different width for the 
annuli. Within a given series, the annuli are all of the same 
width and not overlapping, but there is annulus overlap be- 
tween the series. Thus, the SB profile tabulations from MG03 
include interdependent data points. To avoid such correla- 
tions during our model fitting, we simply adopted only the 
(V+B/I) surface brightnesses derived by MG03 in their nar- 
rowest annuli (AR = 1 "5) for the inner regions of every clus- 
ter, and switched over to take only the widest-annulus results 
(AR = 4") in the outer parts. We discarded altogether the SB 
points estimated from the intermediate annulus widths. Addi- 
tionally, to each annulus that we kept we assigned an average 
(median) radius determined by the local slope of the surface- 
brightness distribution, as described in iKing l 111988) (rather 
than associating each SB point with the straight mean of the 
inner and outer radii of its annulus, as MG03 do). 



At this point we were able to match the ground-based star- 
counts to the combined HST SB profiles, on the V-band mag- 
nitude scale defined by MG03. Background-corrected stel- 
lar number densities N and uncertainties are tabulated for 32 
LMC clusters in the Mackey & Gilmore sample, and for all 
10 of the SMC clusters, in the papers listed in Tabled We 
first converted all of these to uncalibrated "surface bright- 
nesses," -2.5 log N. In many cases, multiple plates with dif- 
ferent exposure times and limiting magnitudes were used by 
the original authors to derive independent estimates of the 
density throughout overlapping ranges of radius R within a 
single cluster. We shifted such partial N(R) profiles by con- 
stants chosen to make the median difference of the various 
log N(R), taken over all radii where any count sets overlap, 
vanish. The single, ground-based profile in -2.5 log N(R) 
that resulted always extended inward to at least R = 100" in 
each cluster, thus overlapping with the MG03 HST /iy data. 
The difference [/xy(/?,)+2.51og N(Rj)] averaged over the radii 
where ground- and space-based data overlap then defined a 
constant C which, when added to -2.5 log A^OR), brings the 
older starcount data onto the V-band surface brightness scale 
of Mackey & Gilmore. In determining these "calibrations," 
we never included the ground-based density at the innermost 
or outermost radius of any count set, as these are the radii po- 
tentially most susceptible, respectively, to crowding and back- 
ground errors. 

With these combined, parsed, and e xtended uv(R) pro- 
files in han d, and in effect calibrated by iMackev & G ilmore 
(2003a be), we found it necessary to re-examine this 
calibration — the zeropoint of their conversion from HST star- 
counts to V-band fluxes — itself. To check it, we used the 
Hv profiles as published in MG03 to calculate the inte- 
grated (enclosed) magnitude, V(< R), as a function of pro- 
jected radius on scales R < 100" in each cluster in the to- 
tal LMC+SMC+Fornax sample. We then compared these in- 
tegrated profiles against appropriate ground-based aperture 
magnitudes. Photometry within apertures R dp < 100" exists 
in the literature for all but two of the clusters under consid- 
eration, and the sources that we have used are listed in Ta- 
ble [3] In principle, we might have simply used true SB pro- 
files from the literature to compare directly against Mackey 
& Gilmore's surface brightness calibration without any inte- 
gration; but such profiles exist for a much smaller fraction of 
the clusters in this sample, so for reasons of homogeneity we 
took the aperture-photometry approach in all cases. 

Table [4] presents the results of this comparison for all 68 
clusters from MG03. Each cluster has six columns in this 
table: first is the cluster name; next, the aperture size, R ap , as- 
sociated with a ground-based magnitude measurement; then 
the V magnitude from the literature and the specific source 
of this V(< R ap ) combination; next, the integrated V magni- 
tude within R dp implied by Mackey & Gilmore's SB profile 
as published; and finally, the difference A/iy of the ground- 
based aperture magnitude minus the MG03 V magnitude. 

For the majority of clusters, the V-band magnitudes implied 
by the MG03 surface-brightness calibration are significantly 
fainter — sometimes by a full magnitude or more — than any 
independent aperture photometry. This appears to be related 
to the fact that, in performing their number counts, MG03 
necessarily had to mask out some bright stars, whether be- 
cause of issues with saturation or scattered light, or because 
they disturbed the smoother overall distribution of the under- 
lying, more numerous fainter stars in the cluster. The flux 
from these stars was not added back in when the stellar counts 
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TABLE 1 

Published Starcount Data for Magellanic Cloud and 
Fornax Clusters 



Source 



Clusters 



LMC: 

Kontizas, Chrvsovereis, & Kontizas (1987) 
Kontizas, Hadiidimitriou, & Kontizas (1987) 
^^^^^^g^^^^^^^^^&mtizas Tim) 

Elson. Fall. & Freeman (1987) 

Mackev & Gilmore (2003a) 



SMC: 

Kontizas, Danezis, & Kontezis (1982) 
Kontizas & Kontizas* 7T%83) 
Mackev & Gilmore (2003b) 

Fornax Dwarf Spheroidal: 
Mackev & Gilmore (2003c) 



NGC 1711, 1786, 1835, 1847, 1850, 1856, 2019, 2100 
NGC 1777, 1868; Hodge 14; SL 842 
NGC 1754, 1805, 1898, 2031, 2121, 2136, 2173, 2210, 
2213,2231 

NGC 1818, 1831, 1866, 2004, 2156, 2157, 2159, 2164, 

2172, 2214 
all of the above, plus: 

NGC 1466, 1651, 1718, 1841, 1860, 1916, 1984, 2005, 
201 1, 2153, 2155, 2162, 2193, 2209, 2249, 2257; 
Hodge 4, 11; R136 (30 Dor); SL 663, 855 



Kron 3; NGC 152, 176, 361, 458 
NGC 121,330, 339,411,416 
all of the above 



FORNAX 1,2, 3,4,5 



TABLE 2 

Instrumental Colors of Magellanic Cloud and Fornax Clusters 



Cluster 


Color 


Value 


Cluster 


Color 


Value 


(1) 


(2) 




(3) 


(1) 


(2) 




(3) 


LMC-HODGE11 


(B- 


V) 


0.35 ±0.09 


LMC-NGC2156 


(V- 


-1) 


0.04 ±0.14 


LMC-HODGE14 


(B- 


V) 


0.41 ±0.26 


LMC-NGC2157 


(B- 


V) 


-0.12 ±0.08 


LMC-HODGE4 


(B- 


V) 


0.42 ±0.11 


LMC-NGC2159 


(V- 


-I) 


0.10±0.11 


LMC-NGC1466 


(V- 


-I) 


0.77 ±0.06 


LMC-NGC2162 


(B- 


V) 


0.46 ±0.11 


LMC-NGC1651 


(B- 


V) 


0.49 ±0.11 


LMC-NGC2164 


(V- 


-I) 


0.13 ±0.06 


LMC-NGC1711 


(V- 


-1) 


0.21 ±0.06 


LMC-NGC2172 


(V- 


-I) 


0.12 ±0.12 


LMC-NGC1718 


(B- 


V) 


0.50±0.10 


LMC-NGC2173 


(B- 


V) 


0.49 ±0.10 


LMC-NGC1754 


(V- 


-1) 


1.02 ±0.07 


LMC-NGC2193 


(B- 


V) 


0.50 ±0.14 


LMC-NGC1777 


(B- 


V) 


0.39±0.10 


LMC-NGC2209 


(B- 


V) 


0.09 ±0.20 


LMC-NGC1786 


(V- 


-I) 


0.87 ±0.05 


LMC-NGC2210 


(V- 


-I) 


0.82 ±0.06 


LMC-NGC 1 805 


(V- 


-I) 


0.08±0.10 


LMC-NGC2213 


(B- 


V) 


0.47 ±0.12 


LMC-NGC 1818 


(V- 


-1) 


0.05 ±0.06 


LMC-NGC2214 


(«- 


V) 


0.00 ±0.09 


LMC-NGC 1831 


(b- 


V) 


0.19±0.05 


LMC-NGC223 1 


(B- 


V) 


0.42 ±0.13 


LMC-NGC 1835 


(V- 


-I) 


1.02 ±0.06 


LMC-NGC2249 


(B- 


V) 


0.27 ± 0.08 


LMC-NGC 1841 


(V- 


-1) 


0.88±0.11 


LMC-NGC2257 


(B- 


V) 


0.37 ±0.10 


LMC-NGC 1847 


(B- 


V) 


-0.05 ±0.11 


LMC-R136 


(V- 


-I) 


0.36 ±0.08 


LMC-NGC 1850 


(B- 


V) 


0.18±0.06 


LMC-SL663 


(B- 


V) 


0.29 ±0.21 


LMC-NGC 1856 


(B- 


V) 


0.26 ±0.04 


LMC-SL842 


(B- 


V) 


0.40±0.19 


LMC-NGC 1860 


(B- 


V) 


-0.07 ±0.14 


LMC-SL855 


(B- 


V) 


0.17 ±0.27 


LMC-NGC 1866 


(V- 


-I) 


0.18±0.05 


SMC-KRON3 


(B- 


V) 


0.42 ± 0.09 


LMC-NGC 1868 


(B- 


V) 


0.29 ±0.07 


SMC-NGC121 


(B- 


V) 


0.49 ±0.06 


LMC-NGC 1898 


(V- 


-I) 


0.86 ±0.07 


SMC-NGC152 


(«- 


V) 


0.39 ±0.11 


LMC-NGC 19 16 


(V- 


-I) 


1.12 ±0.04 


SMC-NGC176 


(B- 


V) 


-0.01 ±0.18 


LMC-NGC 1984 


(V- 


-I) 


0.56±0.12 


SMC-NGC330 


(B- 


V) 


-0.06 ±0.07 


LMC-NGC2004 


(B- 


V) 


-0.13 ±0.11 


SMC-NGC339 


(B- 


V) 


0.44 ±0.11 


LMC-NGC2005 


(V- 


-1) 


0.92 ±0.08 


SMC-NGC361 


(B- 


V) 


0.37 ±0.11 


LMC-NGC2011 


(V- 


-I) 


0.29±0.16 


SMC-NGC411 


(B- 


V) 


0.44 ±0.11 


LMC-NGC2019 


(V- 


-1) 


0.94 ±0.06 


SMC-NGC416 


(B- 


V) 


0.36 ±0.08 


LMC-NGC203 1 


(V- 


-1) 


0.28 ±0.06 


SMC-NGC458 


(B- 


V) 


0.07 ±0.11 


LMC-NGC2100 


(B- 


V) 


0.02 ±0.11 


FORNAX 1 


<y- 


-I) 


0.64 ±0.1 2 


LMC-NGC2121 


(B- 


V) 


0.39±0.10 


FORNAX2 


(V- 


-I) 


0.78 ±0.09 


LMC-NGC2136 


(B- 


V) 


0.17±0.11 


FORNAX3 


(V- 


-I) 


0.89 ±0.08 


LMC-NGC2153 


(B- 


V) 


0.30±0.15 


FORNAX4 


(V- 


-I) 


1.06±0.13 


LMC-NGC2155 


(B- 


V) 


0.51 ±0.13 


FORNAX5 


(V- 


-I) 


0.90 ±0.09 



NOTE. — Cluster colors are those implied by the integrated magnitudes derived from surface-brightness 
profiles out to R 100" as published by Mackey & Gilmore < 2003a b c). Colors are strictly instrumental, in 
that the published surface brightnesses have not been corrected for any zeropoint changes in any bandpass. 
For clusters with instrumental (B-V) colors given here, compare with the true colors listed in Column (4) of 
Tabled 
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TABLE 3 

Published Aperture Photometry for Magellanic 
Cloud and Fornax Clusters 



Source Clusters 



LMC: 

NGC 1835, 2019 
NGC2173 

Mackev & Gilmore (2003a) sample excluding SL 663, 855 
NGC 1835, 1841, 1856, 1916, 2121 
NGC 1818, 1831, 1866, 2004, 2156, 2157, 2159, 2164, 
2172, 2214 

Mackev & Gilmore (2003a) sample excluding SL 663, 855 
SMC: 

van den Bergh (1981) full Mackev & Gilmore (2003b) sample 

Gordon & Kron (1983) full Mackev & Gilmore (200 3b) sample 

Fornax Dwarf Spheroidal: 

Hodge (1965. 1969) : Webbink (1985 ) FORNAX 1 

de VaurouTeurs & AbTesTT^TS) FORNAX 2, 3, 4. 5 

Gordon & Kron (1983) FORNAX 3, 4 



TABLE 4 

Zeropoint Offsets Applied to Mackey & Gilmore V-band Surface Brightness 



Cluster 




Vap 


Ref. 


Vap 




Cluster 




Vap 


Ref. 


Vap 








[sec] 


[lit.] 




[MG03] 






[sec] 


[lit.] 




[MG03] 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


LMC-HODGE11 


30.5 


11.93 


1 


12.158 


-0.228 ± 0.046 


LMC-NGC2156 


27.6 


11.38 


3 


12.670 


-1 


.290 ±0.104 


LMC-HODGE14 


31.0 


13.42 


2 


13.571 


-0.151 ±0.084 


LMC-NGC2157 


30.0 


10.16 


2 


10.920 


-0 


.760 ±0.074 


LMC-HODGE4 


19.0 


13.33 


1 


13.489 


-0.159 ±0.080 


LMC-NGC2159 


36.0 


11.38 


2 


12.569 


-1 


.189 ±0.065 


LMC-NGC1466 


30.0 


11.59 


2 


11.963 


-0.373 ±0.040 


LMC-NGC2162 


31.0 


12.70 


2 


12.814 


-0 


.114 ±0.079 


LMC-NGC1651 


50.0 


12.28 


1 


12.254 


0.026 ± 0.063 


LMC-NGC2164 


30.0 


10.34 


2 


11.603 


-1 


.263 ± 0.040 


LMC-NGC1711 


30.0 


10.11 


2 


11.431 


-1.321 ±0.041 


LMC-NGC2172 


36.0 


11.75 


2 


12.976 


-1 


.226 ± 0.069 


LMC-NGC1718 


31.0 


12.25 


2 


12.219 


0.031 ±0.059 


LMC-NGC2173 


75.0 


11.88 


1 


11.901 


-0 


.021 ±0.059 


LMC-NGC1754 


50.0 


11.57 


1 


11.907 


-0.337 ±0.052 


LMC-NGC2193 


19.0 


13.42 


1 


13.582 


-0 


.162 ±0.089 


LMC-NGC1777 


19.0 


12.80 


1 


12.997 


-0.197 ±0.048 


LMC-NGC2209 


34.0 


13.15 


2 


13.940 


-0 


.790 ±0.147 


LMC-NGC1786 


30.0 


10.88 


2 


11.453 


-0.573 ±0.029 


LMC-NGC2210 


34.0 


10.94 


2 


11.492 


-0 


.552 ±0.033 


LMC-NGC 1 805 


30.0 


10.63 


2 


11.945 


-1.315 ±0.073 


LMC-NGC2213 


31.0 


12.38 


2 


12.466 


-0 


.086 ±0.076 


LMC-NGC 1818 


36.0 


9.70 


2 


11.181 


-1.481 ±0.047 


LMC-NGC2214 


30.0 


10.93 


2 


11.408 


-0 


.478 ±0.088 


LMC-NGC 1831 


30.0 


11.18 


2 


11.219 


-0.039 ± 0.040 


LMC-NGC223 1 


22.0 


13.20 


2 


13.516 


-0 


316±0.105 


LMC-NGC 1835 


31.0 


10.17 


2 


10.652 


-0.482 ±0.037 


LMC-NGC2249 


75.0 


11.94 


1 


12.035 


-0 


.095 ±0.060 


LMC-NGC 1841 


93.5 


11.43 


1 


12.164 


-0.734 ±0.087 


LMC-NGC2257 


30.5 


12.62 


1 


12.747 


-0 


.127 ±0.066 


LMC-NGC 1847 


36.0 


11.06 


2 


11.512 


-0.452 ± 0.076 


LMC-R136 


30.0 


8.27 


2 


8.985 


-0 


.715 ±0.078 


LMC-NGC 1850 


25.0 


9.57 


1 


9.721 


-0.151 ±0.060 


LMC-SL663 















LMC-NGC 1856 


31.0 


10.06 


2 


10.042 


0.018 ±0.032 


LMC-SL842 


19.0 


14.15 


1 


14.196 


-0 


.046 ±0.160 


LMC-NGC 1860 


36.0 


11.04 


2 


12.852 


-1.812±0.134 


LMC-SL855 















LMC-NGC 1866 


36.0 


9.73 


2 


10.838 


-1.108 ±0.034 


SMC-KRON3 


31.0 


12.05 


2 


12.198 


-0 


.148 ±0.046 


LMC-NGC 1868 


31.0 


11.57 


2 


11.633 


-0.063 ± 0.044 


SMC-NGC121 


31.0 


11.24 


2 


11.548 


-0 


.308 ±0.034 


LMC-NGC 1898 


20.0 


11.86 


1 


12.463 


-0.603 ±0.049 


SMC-NGC152 


31.0 


12.92 


2 


12.946 


-0 


.026 ± 0.077 


LMC-NGC1916 


22.0 


10.38 


2 


10.868 


-0.488 ± 0.028 


SMC-NGC176 


31.0 


12.70 


2 


13.441 


-0 


.741 ±0.123 


LMC-NGC 1984 


25.0 


9.99 


1 


12.922 


-2.932 ±0.086 


SMC-NGC330 


31.0 


9.60 


2 


10.621 


-1 


.021 ±0.065 


LMC-NGC2004 


36.0 


9.60 


2 


10.380 


-0.780 ±0.090 


SMC-NGC339 


31.0 


12.84 


2 


12.873 


-0 


.033 ±0.079 


LMC-NGC2005 


12.5 


11.57 


2 


12.052 


-0.482 ±0.058 


SMC-NGC361 


31.0 


12.12 


4 


12.690 


-0 


.570 ±0.085 


LMC-NGC2011 


20.0 


10.58 


1 


13.353 


-2.773 ±0.175 


SMC-NGC411 


31.0 


12.21 


2 


12.234 


-0 


.024 ±0.088 


LMC-NGC2019 


36.0 


10.86 


2 


11.314 


-0.454 ± 0.043 


SMC-NGC416 


31.0 


11.42 


2 


11.609 


-0 


.189 ±0.038 


LMC-NGC203 1 


36.0 


10.83 


2 


11.560 


-0.730 ±0.039 


SMC-NGC458 


31.0 


11.73 


2 


11.832 


-0 


.102 ±0.092 


LMC-NGC2100 


30.0 


9.60 


2 


10.184 


-0.584 ±0.101 


FORNAX 1 


70.0 


15.57 


5 


15.395 







LMC-NGC2121 


31.0 


12.37 


2 


12.581 


-0.211 ±0.058 


FORNAX2 


32.5 


13.73 


6 


14.017 


-0 


.287 ±0.050 


LMC-NGC2136 


30.0 


10.54 


2 


10.733 


-0.193 ±0.089 


FORNAX3 


32.5 


12.74 


6 


13.177 


-0 


.437 ±0.047 


LMC-NGC2153 


50.0 


13.05 


1 


13.472 


-0.422 ±0.103 


FORNAX4 


32.5 


13.49 


6 


14.069 


-0 


.579 ±0.090 


LMC-NGC2155 


31.0 


12.60 


2 


12.542 


0.058 ±0.091 


FORNAX5 


32.5 


13.55 


6 


14.009 


-0 


.459 ±0.055 



REFERENCES. — (l)-IBica et alj 17993) : (2)-[ran den Bereh 1981); (3)— Elson. Fall. & Freeman i 1987); (4)— Gordon & Kron i 1983); (5)— Webbink i 198^); (6)— 
Ide Vaucouieurs & Ablej|197d) . 

NOTE. — Six columns for each of 68 star clusters in the LMC, SMC, and Fornax studied by Mackev & Gilmore > 2003a h c, collectively MG03): 
Column (1) — Cluster name. 

Column (2) — Aperture radius, in arcsec, of a ground-based magnitude measurement from previous literature. 
Column (3) — Ground-based V magnitude of the cluster within the radius in Column (2). 
Column (4) — Source of the cited ground-based aperture magnitude. 

Column (5) — V-band magnitude within obtained by integrating the published surface-brightness profile of MG03. 

Column (6) — Zeropoint correction to be added to MG03 V-band surface brightnesses to bring their cluster magnitudes into agreement with the earlier, ground-based data. 



Bernard & Bisav (1974) 
Bernard (1975) 
van den Bergh (1981) 
Gordon & Kron (1983) 
Elson. Fall.TFreeman (1987) 

Bicaetal. (1996) 
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A/^ v = (literature V ap )-MG03 



FIG. 1. — Histogram of zeropoint corrections A/iy to the V-band 
surface bright ness scale defined for LMC, SMC, and Fornax clusters in 
IMackev & Gilmore 1 2003a b c). Values of Afiy are given for individual clus- 
ters in TableRI 



were converted to mag arcsec -2 , so that the V-band surface 
brightnesses published in MG 03 are sys tematically too faint 
(see the discussion in M ackev & Gilmorell2003al especially; 
also, Mackey 2003, private communication). 

Figure^is a histogram of the zeropoint corrections required 
to bring the MG03 surface brightnesses into agreement with 
the independent aperture photometry we have taken from the 
literature. Some of these offsets are rather surprisingly large, 
particularly in a few of the LMC clusters, and the possibility 
exists in principle that the older aperture magnitudes could be 
the ones in error — if, say, a bright foreground star were unwit- 
tingly included in the ground-based aperture but properly ex- 
cluded from the HST analysis. However, a detailed examina- 
tion of several specific cases has shown that the ground-based 
numbers are not the ones at fault (Mackey 2003, private com- 
munication). More general support for the same conclusion 
comes from the fact that, as we describe just below, there is 
typically good agreement between the integrated V(< R) pro- 
file shapes from the MG03 and ground-based data in many 
clusters where multiple aperture magnitudes can be found in 
the literature — suggesting that if the ground-based data were 
erroneously brightened by foreground stars, these would have 
had to be projected essentially onto the very centers of the 
clusters in an inordinate number of cases. And the instrumen- 
tal colors (Table|2j of those clusters for which MG03 publish 
both B and V surface-brightness profiles do not agree with the 
aperture -photometry colors from the literature [Column (4) of 
Table [8] below], even when the V-band zeropoint correction 
that we infer is negligible (implying that the missing V flux 
in MG03 is not generally consistent with a regular population 
of foreground stars). Thus, when fitting models and deriv- 
ing physi cal p arameters for the Magellanic Cloud and Fornax 
clusters ( 34.21 . we always add the (generally negative) offsets 
AjUy listed in Table0]to the V-band profile numbers published 
in MG03. 



Note that for two of the LMC clusters (SL-663 and SL-855) 
we could find no independent V photometry against which 
to compare the MG03 calibration; while for Globular Clus- 
ter 1 in Fornax , vario us es timates of the total V magnitude 
from IWebbinkl fi985l) and iHodg d (U%5l U%9i) differ from 
each other by at least as much as they do from the implied 
MG03 magnitude. In these three cases, we do not apply any 
zeropoint correction to the MG03 surface brightnesses, but 
we recognize that the adopted numbers could be in error by 
tenths of a magnitude in each case. 

For 32 of the clusters, ground-based V magnitudes are 
available in the literature for more than one aperture — 
allowing for further checks in these cases, both of the consis- 
tency between the various old magnitude estimates and of the 
overall shape of the new Mackey & Gilmore profiles. These 
checks are presented in Figure|2] Each panel here corresponds 
to a different cluster, and all points plotted are aperture magni- 
tudes from the sources list ed in Table|3] Measurements of V ap 
from Ivan den Berg hl lll981l) and lBica et alJ Jl996ft — the most 
comprehensive listings — are set apart as a large open square 
in every panel, while the smaller filled circles refer to reliable 
data from any other source. The large open circles, which 
are particularly evident in the Fornax globular clusters, de- 
note published magnitudes which the original authors have 
indicated are especially uncertain. In all panels of Fig. [2] the 
solid line is the integrated magnitude profile derived from the 
V-band surface brightness data of MG03, after a shift bright- 
ward by the offset A/iy from Tabled which is listed in each 
panel. As was mentioned above, these comparisons show gen- 
erally rather good agreement between the shapes of the cluster 
profiles as measured by MG03, and what can be inferred from 
the ground. In the few cases where the two do appear to differ 
at some significant level (e.g., NGC 1916 in the LMC), we 
simply view the newer HST data as providing an improved 
estimate of the true relative density distribution inside these 
distant clusters. 

Table |5] contains the final surface-brightness profiles of the 
68 LMC, SMC, and Fornax clusters to which we fit models 
in l^below. This table is published in its entirety in the elec- 
tronic version of ApJS; only a short excerpt is presented here. 
Each cluster is given several tens of lines in the table. The 
first column of each line is the cluster name, followed by: a 
radius in arcsec (the median o f an annulu s, determined, as 
mentioned above, according to lKingl[l988l) : the logarithmic 
radius; the V-band surface brightness measured at that radius 
after adjusting the magnitude scale of Mackey & Gilmore by 
the zeropoint offset from Tabled the faint and bright limits 
on nv (reflecting the possibility of asymmetric errorbars, es- 
pecially at faint intensity levels); the bandpass from which 
the datapoint originally came (with G denoting ground-based 
number densities scaled as described above); an estimate of 
the V-band extinction towards the cluster (assumed spatially 
constant in each cluster), and its uncertainty; and a flag indi- 
cating whether or not the point was included explicitly when 
fitting the models of 

To determine whether or not any given point should con- 
tribute to the weighting of the model fits — i.e., in setting the 
"fit flag" of Table [5] to 1 or — we first used all the tabu- 
lated V-B//-ground-based SB values for a cluster to produce 
a smoothed (nonparametric) density profile. Any individual 
points that fell more than 2a away from this smoothed ap- 
proximation were assigned a fit flag of 0. 

The V-band extinctions in Table [5] require further expla- 
nation, since direct estimates of the reddening of individual 
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FIG. 2. — Aperture magnitude vs. radius for LMC, SMC, and Fornax clusters from the sample of Mackev & Gilmore 1 2003a b c). Only clusters with more 
than one ground-based V ap found in the literature are shown. Solid line in each panel is the integrated magnitude V(< R) as a function of radius from the SB 
profile published by MG03, after applying the zeropoint correction given in TablePI Dotted lines demark the uncertainties derived from the MG03 profiles. Open 
squares in all panels denote aperture magnitudes taken from either Bica et al. 1 1996) or van den Berah 1 1981), filled circles are aperture magnitudes from any 
other source (see Tablel3*l, and open circles are measurements cited as uncertain by the original authors. 
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TABLE 5 

Calibrated V-band Surface-Brightness Profiles of LMC, SMC, and Fornax Clusters 



Cluster 


R 






(max) 


(min) 


band b 


A v c 


uncertainty 


fit flag d 




[sec] 




[ma 


g arcsec" 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


LMC-HODGE 1 1 


1.06 


0.025 


18.82 


19.21 


18.43 


V 


0.143 


0.077 




LMC-HODGE11 


1.06 


0.025 


19.00 


19.35 


18.65 


B 


0.143 


0.077 




LMC-HODGE 11 


2.37 


0.374 


19.53 


19.78 


19.28 


B 


0.143 


0.077 




SMC-KRON3 


1.06 


0.027 


19.76 


20.10 


19.42 


V 


0.050 


0.081 




SMC-KRON3 


1.06 


0.027 


19.75 


20.07 


19.43 


B 


0.050 


0.081 




SMC-KRON3 


2.37 


0.375 


20.08 


20.26 


19.90 


B 


0.050 


0.081 




FORNAX 1 


1.41 


0.150 


22.84 


23.29 


22.39 


V 


0.220 


0.077 




FORNAX 1 


1.41 


0.150 


22.61 


23.08 


22.14 


I 


0.220 


0.077 




FORNAX 1 


3.16 


0.500 


23.41 


23.60 


23.22 


I 


0.220 


0.077 





NOTE. — A machine-readable version of the complete Tablel^lis published in the electronic edition of the Astrophys- 
ical Journal Supplement Series; only a portion is shown here, for guidance regarding its form and content. 



^Surface brightnesses fly in Column (4) are on our calibrated V-band scale, with the zeropoint corrections to the 
IMackev & Gilmore ■ 2003a b c) numbers already applied (see Tablel4l. Maximum and minimum surface brightnesses 
in Columns (5) and (6) define the errorbars on fly, which are asymmetric in general. The uncertainty in the zeropoint 
correction A/xy from Tablepts a source of systematic error; it is not accounted for in Columns (5) and (6) but is added in 
quadrature to the formal, x*-based uncertainties in all fitted and derived cluster parameters relating to surface brightness. 

^V, B, or / denotes the original bandpass of a Mackev & Gilmore i 2003a b c) surface-brightness datapoint, with B or 
/ data shifted in zeropoint as in Tablel2l and then again by an amount A fly from Tablel^l to match onto our calibrated 
fly scale. G denotes a surface density taken from ground-based starcounts in the literature (TablefTl. transformed to a 
surface-brightness and zeropointed to the fj,y scale of column (4). 

c V-ba nd extinction Av = 3.1 EiB — V), with E(B — V) taken from literature values compiled by Mackev & Gilmore 
12003d) for the Fornax globular clusters but inferred from our population-synthesis modeling in $J3](Table[$] for LMC 
and SMC clusters. Ay must be subtracted from columns (4), (5), and (6); it is spatially constant in each cluster but is 
listed at every radius as a convenience. 

^ 1 = point used in fitting of structural models in ^4] = point not used in model fitting. 



Magellanic Cloud clusters do not exist for the majority in this 
sample. However, in order to accurately determine cluster 
physical parameters — true central surface brightnesses, total 
luminosities and masses, etc. — some knowledge of the extinc- 
tion is clearly required. MG03 assign a single, average ex- 
tinction to all LMC clusters, and another average to all SMC 
clusters; but the reddening is known to vary across these sys- 
tems. We have instead used population-synthesis modeling, 
given an age and a metallicity for each cluster as compiled 
from the literature by MG03, to predict an intrinsic (B-V)o 
color for every object here. Observed (B-V) colors from 
the aperture-photometry literature listed in Table|3]then imply 
reddenings E(B-V), and extinctions Ay = 3.1E(B-V) fol- 
low. More details of this procedure are given in |5] where we 
also use population-synthesis models to produce estimates of 
cluster mass-to-light ratios. In the next subsection, we first 
briefly describe some point s related to our handling of the 
Trager. King. & Diorgovski ( 1995) catalogue data for Milky 
Way globular clusters. 

2.2. Milky Way Globular Clusters 

As was mentioned above, Trager . King:. & D iorgovski 
i ll 9951) have already taken steps, similar to those we have out- 
lined above and applied to LMC, SMC, and Fornax clusters, 
to combine heterogeneous data into unique, calibrated V-band 
SB profiles for 124 Galactic globular clusters. However, their 
results as published do not include any estimates of uncer- 
tainty in the individual surface-brightness points. The pro- 
files defined above for the MG03 cluster sample do include 
errorbars, and we use these during our model fitting in §4 to 
estimate the uncertainties in all basic and derived physical pa- 
rameters of the clusters, via a standard analysis of x 2 varia- 
tions over grids of model fits. We would like to proceed in 
the same way with the Milky Way GC sample, and thus we 



ha ve attempted to estimate errorbars for the /jy values given 

bv lTrager. King. & Diorgovski i ll 9951). 

In lieu of absolute errorbars. iTrager. King. & Diorgovskil 
give each point in their brightness profiles a relative weight 
w< £ [0, 1]. As they discuss, these weights were assigned 
"by eye" to reflect the authors' judgement of the overall 
quality of the source dataset. They are therefore not con- 
nected rigorously to relative errors, and their precise mean- 
ing is left somewhat open to interpretation. Initially, we pro- 
ceeded under the natural assumption that the weights were 
proportional to the inverse square of the surface-brightness 
uncertainties, or u, oc 1 / yfw~i. When fitting models by 
minimizing an error-weighted \ 2 statistic, this led to some 
cases where the best-fit model parameters were unduly influ- 
enced by just a few d iscrepant points with low weights from 
Trager, King, & Diorgovski. Thus, we decided instead to de- 
weight such points even further by adopting the heuristic pre- 
scription that the relative surface-brightness errorbars grow as 
07 oc 1 /wj. 

Give n this choice, we e stimated the un certainties as fol- 
lows. Trager, King, & Diorgovski ( 1995) also tabulate the 
value at each radius in each cluster, of interpolating (Cheby- 
shev) polynomials that provide reasonably accurate, model- 
independent approximations to the overall SB profile. We 
assume that the reduced \ 2 P er degree of freedom for these 
polynomial fits is e xactly 1 for every cluster. T ypically, the 
polynomials fit by Trager. King. & Diorgovski are of third 
order, so for a surface-brightness profile with N datapoints 
there are A^-4 degrees of freedom. Then, with our pre- 
scription er,' = tT M /w,-, where ov is a different constant for 
each GC, and writing C{R/) for the value of the polyno- 
mial fit at each radius in any one cluster, we have that 
jhT.U^/^lnCiRd-MRi)] 2 = 1. We have solved this 
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TABLE 6 

Base Errorbars Assigned to Milky Way Globular-Cluster Surface 
Brightnesses in Catalogue of Tracer et al. (1995) 



Cluster 




Cln stpr 




Cln stpr 




ClllQtfM" 
V^lLla LCI 


°> 


AMI 


0.300 


NGC5986 


0.081 


NGC6397 


0.199 


NGC6752 


0.173 


ARP2 


0.200 


NGC6093 


0.171 


NGC6401 


0.119 


NGC6760 


0.207 


HP 


0.248 


NGC6101 


0.182 


NGC6402 


0.137 


NGC6779 


0.207 


IC1276 


0.259 


NGC6121 


0.180 


NGC6426 


0.238 


NGC6809 


0.197 


IC4499 


0.218 


NGC6139 


0.136 


NGC6440 


0.110 


NGC6864 


0.138 


NGC104 


0.082 


NGC6144 


0.155 


NGC6441 


0.076 


NGC6934 


0.173 


NGC1261 


0.187 


NGC6171 


0.143 


NGC6453 


0.150 


NGC6981 


0.223 


NGC1851 


0.249 


NGC6205 


0.128 


NGC6496 


0.307 


NGC7006 


0.165 


NGC1904 


0.094 


NGC6218 


0.182 


NGC6517 


0.098 


NGC7078 


0.178 


NGC2298 


0.147 


NGC6229 


0.151 


NGC6522 


0.140 


NGC7089 


0.107 


NGC2419 


0.120 


NGC6235 


0.188 


NGC6528 


0.283 


NGC7099 


0.115 


NGC2808 


0.124 


NGC6254 


0.137 


NGC6535 


0.348 


NGC7492 


0.253 


NGC288 


0.158 


NGC6256 


0.197 


NGC6539 


0.235 


PALI 


0.202 


NGC3201 


0.206 


NGC6266 


0.126 


NGC6541 


0.201 


PAL 10 


0.075 


NGC362 


0.097 


NGC6273 


0.124 


NGC6544 


0.205 


PAL 11 


0.220 


NGC4147 


0.242 


NGC6284 


0.122 


NGC6553 


0.162 


PAL 12 


0.395 


NGC4372 


0.428 


NGC6287 


0.158 


NGC6558 


0.117 


PAL 13 


0.454 


NGC4590 


0.176 


NGC6293 


0.146 


NGC6569 


0.240 


PALM 


0.140 


NGC5024 


0.153 


NGC6304 


0.148 


NGC6584 


0.148 


PAL2 


0.359 


NGC5053 


0.244 


NGC6316 


0.139 


NGC6624 


0.130 


PAL3 


0.120 


NGC5139 


0.142 


NGC6325 


0.162 


NGC6626 


0.137 


PAL4 


0.129 


NGC5272 


0.191 


NGC6333 


0.169 


NGC6637 


0.093 


PAL5 


0.190 


NGC5286 


0.128 


NGC6341 


0.119 


NGC6638 


0.111 


PAL6 


0.098 


NGC5466 


0.120 


NGC6342 


0.165 


NGC6642 


0.140 


PAL8 


0.220 


NGC5634 


0.220 


NGC6352 


0.269 


NGC6652 


0.190 


TERZAN1 


0.142 


NGC5694 


0.120 


NGC6355 


0.152 


NGC6656 


0.173 


TERZAN2 


0.157 


NGC5824 


0.106 


NGC6356 


0.074 


NGC6681 


0.186 


TERZAN5 


0.129 


NGC5897 


0.118 


NGC6362 


0.136 


NGC6712 


0.218 


TERZAN6 


0.099 


NGC5904 


0.104 


NGC6366 


0.217 


NGC6715 


0.117 


TERZAN7 


0.236 


NGC5927 


0.089 


NGC6380 


0.173 


NGC6717 


0.215 


TERZAN9 


0.128 


NGC5946 


0.136 


NGC6388 


0.101 


NGC6723 


0.212 


TON2 


0.162 



NOT E. — a M is a constant for each Galactic globular cluster in the catalogue of Trager. King. & Dioraovski 
1199 ."jit , used to estimate errorbars on individual surface-brightness datapoints. Given a relative weight Wi (0 < 
vv, < 1) from Trager et al. for each datapoint in a cluster, we define the uncertainty in [iy to be cr, = cr^/vv,. 



identity for the "base" errorba r tr„ of each of the 124 cl usters 
in the collection of Trag er. King. & Diorgovskil ( 119951) . and 
reported the results in Table [S] Again, the uncertainty in /xy 
at any single radius in any one cluster is taken as cr^/vv,, 
assumed to be symmetric about the measured surface bright- 
ness. 

For two globular clusters — Palomar 10 and Terzan 7 — 
Trager et al. only present uncalibrated surface brightnesses, 
jity profiles relative to an unknown central value. Subse- 
quently, the central brightnesses of these objects have been 
determined, and they are tabulated in the catalogue of Har- 
ris (1996): /j,yfi = 22.12 for Palomar 10, and /iyo = 20.69 for 
Terzan 7. We have used these values as zeropoint "correc- 
tions" to the Trager et al. data, and treated them subsequently 
in the same fashion as our zeropoint shifts to all of the Mackey 
& Gilmore LMC/SMC/Fornax cluster surface brightnesses. 

Aside from these two small points, we have proceeded with 
modeling the Trager et al. data as published, with no further 
embellishments nor any attempted updates. There are another 
26 globular clusters included in the more recent catalogue of 
lHarrisI (119961) . but we have not made any attempt to model 
their density structures: the raw profiles of these additional 
clusters have not been gathered into a uniform collection on 
par with the Trager et al. database, and indeed many of them 
are rather obscure and not studied well enough for our pur- 
poses in the first place. For completeness, we list them in Ta- 
ble0as objects that we have not fit with structural models in 



|4] We do, however, include most of them in the population- 
synthesis modeling of |3j since only estimates of their metal- 
licities are required there. 

We originally attempted to fit all 124 globulars from 
iTrager. King. & Diorgovskil (119951) with the variety of mod- 
els that we also fit to the massive clusters in the Magellanic 
Clouds and Fornax. In doing so we found that for several 
clusters the effective (projected half-light) radius R/, returned 
by the fits was of order or larger than the projected radius 
7?i ast of the outermost SB datapoint tabulated by Trager et 
al. In particular, there are 34 GCs for which Ri,/Ri as t > 0.9, 
signifying either that the observations are too sparse to sig- 
nificantly constrain important aspects of the models, or that 
the models we fit are simply poor descriptions of the data. 
In fact, these 34 objects include a majority of those desig- 
nat ed as core-collapse or "possible" co re-collapse candidates 
bv ITrager. King. &D iorgovski (1995), a nd subsequen tly in 
the catalogues of Diorgovski ( 1993) and Harris ( 199(1 ). As 
Trager et al. make very clear, their own estimation of iKingl 
( 1966) model parameters for such globulars are essentially 
rough guesses and not quantitatively trustworthy. Our fits to 
the 34 GCs with Rh/R\- dst > 0.9 are likewise rather uncertain, 
and we do not present any detailed results for them in |4]and 
later sections. Again, however, we name the clusters in Table 
0for completeness, and we include them in the population- 
synthesis analysis that we describe next. 3 After making this 

3 In i|4.3.2l it will be seen that there are also 4 LMC clusters in the MG03 
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TABLE 7 

Milky Way Globular Clusters not fit by Structural Models 



Reason 



Clusters 



Clusters in Harris (1996) but 

not in^rag^^^mg i & Diorgovski (1995) 



1636-283; 2MS-GC01, GC02; AM4; BH176; DJORG 1, 2; 
E 3; ERIDANUS; ESO-SC06; IC 1257; LILLER 1; LYNGA 7; 
NGC 4833, 6540, 6749, 6838; PAL 15; PYXIS; RUP 106; 
TERZAN 3, 4, 8, 10, 12; UKS 1 



Clusters in Trager, King, & Diorgovski (1995) 
with large fitted R;, (see text) 



NGC 4372, 5927, 5946, 6144, 6256, 6284, 6293, 6304, 
6325, 6342, 6352, 6355, 6380, 6401, 6426, 6453, 
6517, 6522, 6544, 6558, 6624, 6626, 6642, 6717; 

HP; PAL 6, 8, 13; TERZAN 1, 2, 5, 6, 9; TON 2 



Core-collapsed clusters 



NGC 6397, 6681, 6752, 7078, 7099 



cut, five other core-collapsed globulars remain in the Trager 
et al. sample. These are known a priori not to be properly 
described by King models (or others with constant-density 
cores) and so we remove them as well from consideration for 
structural modeling, but list them in Tableland include them 
in our population-synthesis modeling. 

3. POPULATION-SYNTHESIS MODELS: REDDENINGS AND 
MASS-TO-LIGHT RATIOS 

Before proceeding with the structural and dynamical mod- 
eling of the clusters in our Magellanic Cloud, Fornax, and 
Milky Way sample, we discuss some aspects of population- 
synthesis modeling. We have folded this into our derivation 
of physical cluster parameters from model fits to the cluster 
SB distributions. There are two main issues at hand that ne- 
cessitate these considerations: 

First, in order to move from a description of the observed 
surface brightness profile of a cluster to one of the true lu- 
minosity density profile, we need information on the extinc- 
tion towards the cluster. For each of the Milky Way globular 
clusters that we model, an estimate of the fore ground redden- 
ing E(B-V) is available in the catalogue of lHarrisI |l996); 
the extinction follows directly, and the required corrections 
are straightforward. The same is true of the globulars in For- 
nax: lMackev & Gilmorel(l2003cll . for example, have tabulated 
measurements from the literature of the reddenings E(B-V) 
of each cluster. But, as was mentioned above, there is no 
such information in the literature for many of the 63 LMC 
and SMC clusters being analyzed here. Aperture measure- 
ments (ground-based) of the clusters' (B-V) dp colors do ex- 
ist, however (see the references in Table [3j, as do individual 
estimates of their ag es and me tallicities (compiled from the 
literature, again, by Mackev & Gilmore 2003a b). We there- 
fore use the cluster ages and [Fe/H] values to derive an ex- 
pected intrinsic (B - V)o from population-synthesis models, 
and then E(B-V) = (fi- V) ap -(fi- V) and A v = 3.1E(B-V). 
The fact that direct reddening measurements do exist for the 
Milky Way and Fornax globular clusters allows us to perform 
the same analysis on them, and then compare our "theoreti- 
cal" E(B-V) values to the known ones — a valuable check on 
the method, albeit only in the extreme of old ages. 

Second, to go from a description of the de-reddened lu- 
minosity density profile of a cluster to its mass density dis- 

sample which have Ri,/R[a S t > 1 and particularly uncertain fit extrapolations. 
In these cases, however, this is because the clusters have rather large cores 
rather than any hint of a post-collapse morphology. Given this, and the fact 
that they have not been fully modeled before, we do report all of our results 
for these objects. 



tribution, knowledge of an average mass-to-light ratio is re- 
quired. Given a sample of clusters of a roughly common age 
(and a common stellar IMF), the mass-to-light ratio typically 
can also be treated as roughly constant, at least to within a 
factor of order unity which depends on [Fe/H] differences. 
Luminosity differences are then essentially proportional to 
mass differences, and detailed knowledge of the exact mass- 
to-light ratio is not critical to the accurate discernment of rela- 
tive trends in physical cluster properties (such as mass-radius, 
mass-velocity dispersion, or other fundamental-plane corre- 
lations). In our case, however, we aim ultimately to com- 
pare these sorts of correlations for the old globular clusters 
in the Milky Way and other galaxies (ages r ~ 13 Gyr for 
the most part), against those for the much younger massive 
clusters in the Magellanic Clouds (t < 10 9 yr and as young 
as r ~ 3 x 10 6 yr for R136=30 Doradus). Trends or corre- 
lations in luminosity then reflect a complex mix of age and 
mass effects which must be separated to make sense of the 
physical situation. Ideally, we would have liked to use mea- 
surements of the stellar velocity dispersions in the clusters to 
compute their mass-to-light ratios directly; but, as for the red- 
denings, such measurements exist for only a handful of the 
young LMC/SMC clusters, and they tend to be highly uncer- 
tain. Even in the Milky Way globular cluster system, reliable 
velocity-dispersion measurements and dynamical M/L ratios 
exist for fewer than half of the 85 clusters that we model here. 
Therefore, we also use population-synthesis models to define 
a V-band mass-to-light ratio for every cluster in our total sam- 
ple. In 5j5]we compare these predicted values with the dynam- 
ical mass-to-light ratios that can be computed for the minority 
of (mostly old) clusters with measured velocity dispersions. 

The population-synthesis model that we have chosen to use 
to derive cluster reddenings and mass-to-light ratios is that 
of iBruzual & Charlotl (120031) . using the Padova 1994 stellar- 
ev olution tracks an d assuming a stellar IMF following that 
of Chabrierj J2003I) fo r the Galactic disk population. The re- 
sults are presented in 33.21 However, since they are so central 
to the establishment of all final, physical (mass-based) clus- 
ter properties and to the definition of systematic interdepen- 
dences between these properties, we have also carried through 
the full suite of calculations with the alternate population- 
synthesis code PEGASE dFioc & Rocca-Volmerangel 119971 
version 2.0) and unde r dif ferent assumptions on the form 
of the stellar IMF. In 33.11 then, we first present point-by- 
point comparisons between the E(B-V) and Ty =M/Ly val- 
ues predicted by various combinations of codes and IMFs. 
Note, however, that we always assume that every clus- 
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FIG. 3. — Top panel: Intrinsic color as a function of age for single- 
burst stellar populations of three different metallicities, as predicted by the 
population-synthesis code of Bruzual & Charlol 1 2003 ) assuming the disk- 
star IMF of Chabriei| 120031) . Bottom panel: Difference in (B - V) for 
the same model clusters a s predicte d by the PEG ASE (v.2.0) population- 
synthesis code of Fioc & Rocca-Volmerange 1 1997). 



ter is a single-age population formed instantaneously in 
one coherent burst of star formation, and that all clus- 
ters, of any age or metallicit y in any galaxy, share a com- 
mon stellar IMF. Also, bo th Bruzual & Chariot! (120031) and 
Fioc & Rocca-Volmerange (1997) have coded prescriptions 
for mass loss over time due to stellar-evolution debris (from 
winds and supernovae, essentially) that is assumed to be 
swept out of a cluster. When using either code, we always 
employ these prescriptions. 

3.1. Comparison of Codes and Stellar IMFs 

The IChabrierl 120031) disk-s tar IMF t hat w e adopt for 
our primary calculations is a Salpeter (1955) power law 
(dN I dm oc to -2,35 ) for stellar masses m > 1M , and a much 
flatter, lognormal distribution below m < 1M . Figures [3] 
and 0] show, in their larger (upper) panels, the intrinsic col- 
ors (B^y]r2_andJ£;bmdjnass-to-light ratios Ty predicted by 
the Bruzual & Chariot (2003) code for a single-burst (or "sim- 
ple") stellar population with this IMF, as functions of cluster 
age for three set heavy-element abundances roughly spanning 
the range appropriate for our combined sample of globular 
and young massive clusters. The smaller (bottom) panels in 
these figures show the differences in (B-V)o and log Ty as 
predic ted by the PEGASE code of lFioc & Rocca-Volmerangel 
( 1997) given the same IMF and the same metallicities. Aside 
from the problematic region around ages 10 6 < t < 10 7 yr, 
there is generally good agreement between the two codes. 

Figure |3 expands somewhat on these plots, showing (B — 
V)o and Yy predicte d as func ti ons of cluster [Fe/H] at var- 
ious ages, all for a Chabrier (120031) disk IMF. The loga- 
rithm of cluster age is labeled in the upper panel of each pair 



FIG. 4. — Top panel: V-band mass-to-light ratio as a function of age for 
single-burst stellar populations of three different metallicities, as predicted 
by the population-synthesis code of Bruzual & Chariot 1 2003) assuming the 
disk-star IMF of Chabriei 1 2003). Bottom panel: Difference in log Yy for 
the same model clusters as predicted by the PEGASE (v.2.0) population- 
synthesis code of lFioc & Rocca-Volmerangel 119971) . Both codes allow for 
mass loss from the clusters due to the evacuation of stellar-wind and super- 
nova debris over time. This amounts to a ~ 30% reduction in total cluster 
mass after a Hubble time. 



in this figure. Filled circ les are the quantities predicted by 
Bruzual & Chariot (2003) at the six metallicities for which 
their code has stellar-evolution tracks. The solid curve join- 
ing these points is a spline interpolation, which we use to infer 
(B-V)o and Ty for clusters with metallicities between the few 
that are explicitly calculable (at any age) with the Bruzual- 
Charlot code. The open squares in every panel are the syn- 
thetic colors and mass-to-light ratios predicted by PEGASE 
v.2.0, again at the six metallicities for which that code has ex- 
plicit stellar-evolution tracks (these metallicities, and indeed 
the tracks themselves, are the same as those used by Bruzual 
& Chariot). The dashed lines joining the squares are linear in- 
terpolations for determining quantities at intermediate [Fe/H]. 
(The PEGASE code does this interpolation internally, given 
any arbitrary [Fe/H] specified by the user.) 

More detailed discussion of the comparison between the 
PEGASE and Bruzual-Charlot codes — including the difficul- 
ties at young ages r < 10 7 years, which f ortunately are not 
relevan t to us in general — may be found in Bruzual & Chariot 
(2003). The plots we have presented have vertical-axis scales 
chosen deliberately to emphasize the differences between the 
codes for one IMF. Evidently, although there are some ex- 
ceptions, these differences are generally slight: at the level 
of < 0.05 mag in (B-V)o and < 10% in Ty. Again, we have 
adopted the Bruzual & Chariot code to define the intrinsic col- 
ors and mass-to-light ratios for our clusters, and we can pro- 
ceed with some confidence that this choice is not introducing 
sizeable systematic errors in our subsequent analyses. 
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FIG. 5. — Intrinsic (B — V) colors and V-band mass-to-light ratios as functions of metallicity for single-burst stellar populations at a range of fixed ages. Filled 
circles and solid lines in every panel refer to predictions from the population-synthesis code of Bruzual & Chariot (2003); open squares and dashed lines refer to 
predictions from the code of Fioc & Rocca-Volmeranse 1 1997, PEGASE v.2.0). The logarithm of the cluster age is given in the upper (B-V)q panel of each pair. 
All calculations employ the disk-star IMF o flChabrieiH 2003). 
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FIG. 6. — Top panel: IMF dependence of intrinsic color (B — V)o of 
a single-burst stellar population in the population-synthesis code PEGASE 
(v.2.0) i Fioc & Rocca-Volmeranee 1997). Bottom panel: IMF dependence 
of the logarithm of V-band mass-to-light ratio log Ty in the PEGASE code. 
Results are shown as differences of the computed quantities using the vari- 
ous IM Fs indicated, minus the "reference" values obtained with a Chabrier 
1 2003) disk-star IMF. Bold lines show the differences as functions of cluster 
age for a metal abundance Z = 0.004 = O.2Z0; lighter lines, when shown, 
correspond to Z = 0.0004 and Z = 0.02. 



The next question concerns the implications of our choice 
of the Chabriefl ( 120031) disk-star IMF for the population- 
synthesis calculations. The Bruzual-Charlot code offers only 
the choice between this option or a pure Salo eterl (119551) 
power law at all stellar masses (which is well known by now 
to be an incorrect description of the true IMF below ~ 1 M , 
but is still commonly used for reference calculations). By 
contrast, the PEGASE code of Fioc & Rocca- Volmerange can 
produce models for any user-defined IMF. We have therefore 
used PEGASE to calculate (B-V)o and Ty as functions of age 
for clusters of three different metal abundances (Z = 0.0004, 
Z = 0.004, and Z = 0.02, as in Figs.|3]and0J, given three stellar 
IMFs besides our preferred Chabrier ( 2003 ) disk-star distribu- 
tion. 

The alternate IMFs that we have examined are those of 
Salceterl (119551) ("a single power law, dN /dm oc m~ 235 at all 
masses); Krouna. Tout. & Gilmore ( 1993) (a three-part piece- 
wise power law which, with dN/dm oc m~ 2 J above m > 1 M , 
is significantly steeper than either Salpeter or Chabrier at high 
masses, but significantly shallower than Salpeter — and not 
substantially different from Chabrier's disk IM F — at lower 
masses); and the "globular cluster" IMF of Chabrier (2003) 
(which is nearly the same as his disk-star IMF but for a some- 
what narrower and displaced lognormal peak below 1M ). 
Figure |6] shows the differences in the predicted (B-V)o and 
logarithmic Ty predicted by PEGASE for these IMFs vs. the 
Chabrier (2003) disk-star function. 

The bold solid lines in both panels of Figure [6] compare 



the Chabrier disk and Sal peterl (119551) IMFs: the bold dashed 
lines, the Chabrier disk and Kroupa. Tout. & Gilmore (Tf99l 
IMFs; and the bold dotted lines, the Chabrier disk and "glob- 
ular cluster" IMFs. All of these bold lines correspond to a 
cluster metal abundance Z = 0.004 = 0.2Z Q ; the thinner lines 
around them refer to Z = 0.0004 = 0.02 Z and Z = 0.02 = Z Q 
(for clarity, in the upper panel the lines for different metallici- 
ties are only shown for the Chabrier disk vs. Kroupa et al. IMF 
comparison). Clearly, the choice of IMF has little bearing on 
the intrinsic (B - V)o color of a cluster at a given age, with 
changes of less than a few hundredths of a magnitude — within 
the differences between the PEGASE and Bruzual-Charlot 
codes — typically being implied. Given that the color is the 
ratio of the total cluster luminosity in two bandpasses, and 
the two codes employ essentially identical treatments of stel- 
lar evolution, its robustness against even major changes in the 
IMF is expected. 

The V-band mass-to-light ratio is, however, naturally more 
sensitive to details of the IMF. The Salpeter IMF gives Ty 
consistently higher than the Chabrier disk IMF (by factors of 
~ 60% - 75%), because of the much higher proportion of the 
total cluster mass that the former distribution places in very 
faint, low-mass stars. As we mentioned just above, however, 
it is known that the Salpeter IMF is not a good desc ription of 
the re al distribution at low stellar masses (see, e.g., Chabrier 
2003). The difference in Ty between Chabrier's disk and 
"globular cluster" IMFs is essentially negligible until very old 
ages, where the GC IMF would predict slightly lower mass- 
to-light ratios. This is because the two IMFs are identical 
at stellar masses m > 1M , where both take on the Salpeter 
power-law shape. It is only after these relatively massive stars 
have all evolved significantly that the slight difference be- 
tween the two IMFs at lower stellar masses has a measurable 
effect. Even then the effect is small, however, and in any case 
Chabrier's disk IMF is better constrained empirically than his 
GC IMF — both arguments for our use of the disk-star IMF to 
make predictions for all of our clusters, GCs as well as the 
massive young objects in the LMC and SMC. 

By con t rast, th e difference betwe en Ty as computed for the 
Chabrier (2003) disk IMF vs. the iKroupa. Tout. & Gilmorel 
(1993) IMF is rather dramatic and, more worrisome, varies 
systematically with cluster age. This is a direct result of the 
much steeper slope specified by the Kroupa et al. IMF for 
stellar masses above 1 M . For a given total cluster mass, the 
Kroupa et al. model initially puts a much smaller fraction into 
massive (and bright) stars, yielding significantly higher mass- 
to-light ratios at young ages. Over time, Ty for the Kroupa et 
al. IMF approaches Ty for the Chabrier IMF, because the two 
IMFs are effectively the same at stellar masses m < 1 M — 
which, of course, contribute all the cluster light at ages of 10 
Gyr and more. 

Clearly, the choice between these two IMFs has a di- 
rect impact on the mass-to-light ratios and all dependent 
physical parameters of clusters younger than t < 10 10 yr — 
which is to say, all of the young LMC and SMC clusters 
that we model. In principle, given the strong systemat- 
ics in the bottom panel of Fig. [6] choosing between these 
two particular IMFs might even influence the ultimate in- 
ference of correlations betwe en (mas s -based ) cluster prop- 
erties. We much prefer the Chabrier (2003) disk-star IMF, 
which has stronger and more recent empirical support than 
the older Kroupa et al. function. Bu t it is important to note 
that iMackev & Gilmorel (l2003albld) adopted the Kroupa el 
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al. IMF for their own population-synthesis modeling. Our es- 
timates below of the young LMC/SMC cluster mass-to-light 
ratios therefore differ significantly from those presented by 
Mackey & Gilmore. However, the lower Yy implied by our 
choice of IMF combines with our generally brighter V-band 
surface-brightness scale ( 32. 1> to produce cluster masses and 
mass densities that often differ less strongly from those in 
MG03. 

3.2. Results for our Cluster Sample 

Table [8] (a sample of which can be found at the end of this 
preprint, and which can be downloaded in full from the elec- 
tronic edition of the Astrophysical Journal Supplement Series) 
presents in full the results of our population-synthesis model- 
ing of the 53 LMC clusters, 10 SMC clusters, and 5 Fornax 
globular clusters in the combined MG03 database, as well as 
for 148 Galactic globular clust ers with [Fe/H] values given in 
the catalogue of lHarris| ([T996) (that is, 58 GCs in addition to 
the 90 that we fit with structural models in §3). In every case, 
to predict an intrinsic (B - V)o color and average (global) Ty 
ratio, we require only an estimate of the cluster age and metal- 
licity. In the LMC, SMC, and Fornax cases, these quantities 
have been either derived or recovered from the literature by 
MG03, and we have generally taken their tabulated values, 
with the single exception that we assign an age of r = 13 ±2 
Gyr (logr = 10.11 ±0.07) to every cluster that they list as 
older th an 13 Gyr. In the Milky Way, we have taken [Fe/H] 
from the Harris ( 1996) catalogue, and assigned a single age of 
13 ± 2 Gyr to all clusters. (Although an age spread of a few 
Gyr is known to exist in the Galactic GC system, at such an 
old average age this makes little difference in the population- 
synthesis colors and mass-to-light ratios; see Figs. [3] and 0J) 
The second and third columns of Table [8] list these ages and 
metallicities. Column 4 gives the observed (B — V) color of 
the cluster whenever such a measurement exists (taken either 
from the aperture -photometry sources in Table [5] above, or 
from IHarrisl 1 1 9961 for the Milky Way globulars)). The rest 
of the two lines for each cluster list the intrinsic colors and 
M/Ly ratios obtained from each of six combinations of two 
pop ulati on-synthesis codes and four stellar IMFs discussed 
in 33.11 I t is column (6 ) of Table HI — the com bination of 
Bruzual & Chariot! (120031) model with Tchabrierl (120031) disk- 
star IMF — that we draw on for all cluster modeling that fol- 
lows. All uncertainties on the population-synthesis quantities 
in Table|H]follow directly from uncertainties in the cluster ages 
and [Fe/H] values (for which errorbars are cited by MG03 or 
IHarrisl 19961 but not reproduced here). 

We are not aware of any previous, comprehensive calcula- 
tion of theoretical mass-to-light ratios for Galactic globular 
clusters. In Fig. Q we show the distribution of population- 
synthesis Ty values for the Milky Way GC system, taken 
from Column (6) of Table [8] Note the strong concentration 
of clusters at T v op ~ 1 .9M Q L^ v , which is to be compared to 

the average dynamically determined (T y yn ) = 1 .45 M^L^y 
for a much smaller sample of clusters (McLauahlin 2000). 
The tail towards higher T v op values in Fig.0is a direct reflec- 
tion of the metal-rich (bulge-cluster) tail in the distribution of 
GC metallicities in the Galaxy. A more detailed, cluster-by- 
cluster comparison of population-synthesis and dynamically- 
measured mass-to-light ratios is given in |5]below, following 
the bulk of our model fitting in 50] 

In the top panel of Fig. [8] we compare our adopted 
population-synthesis mass-to-light ratios for the 




FIG. 7. — V-band mass-to-light ratios of 148 Galactic globular clusters 
as predicted by the population-synthesis code of Bruzual & Chariot 1 2003) 
using the disk-star IMF of Chabrier 1 2003). Values for individual objects are 
in Column (6) of Table 151 A common age of 13 ± 2 Gyr has been assumed 
for all clusters, so that the spread in Yy° p directly reflects the metallicity 
distribution of Galactic GCs given [Fe/H] values taken from Harris 1 1996). 
The horizontal errorbar shows the rms uncertainty in the predicted Ty . 



LMC+SM C+Fornax cluster sample, with those com- 
puted by Mac kev & Gilmorel d2003alblcl) . Ours are always 
lower than theirs, in part because of the much s teep er stellar 
IMF that MG03 assumed — as discussed in 33. H — and in 
part because they appear not to have allowed for the loss of 
cluster mass through massive-star ejecta in their application 
of the PEGASE model (a ~ 30% cumulative effect over 
a Hubble time, in either the Bruzual-Charlot or Fioc & 
Rocca-Volmerange codes). Further still, MG03 give the ages 
of some of their oldest LMC clusters — the genuine globulars 
there — as greater than 13 Gyr, but we have fixed such ages 
at t = 13 ± 2 Gyr, resulting again in lower T v ° p ratios. We 
note in passing that the leftmost point in this plot — the cluster 
with the lowest mass-to-light ratio in MG03, and the largest 
errorbar in our own assessment — is Rl 36=30 Doradus, which 
at t ~ 3 x 10 6 yr is the youngest cluster in our entire sample 
and the one most susceptible to uncertainties in the current 
population-synthesis models. 

Finally, the bottom panel of Fig.Elplots the difference be- 
tween our population-synthesis derived E(B - V) = [(B — V) - 
(B-V)o] — that is, the difference of Columns (4) and (6) in 
Table [8] — and the published reddenings for individual clus- 
ters in the LMC, SMC, Fornax, and the Milky Way. This 
plot is dominated by the Galac tic GC sample , where the pub- 
lished E(B-V) are those from lHarria 111996). These are gen- 
erally in very good agreement with our model calculations; 
the strongest outliers are Palomar 6 and Terzan 5, which have 
large E(B-V) = 1.46 and 2.15, respectively. This reflects well 
both on the calibration of the population-synthesis models 
against old, single-burst stellar pop ulations and o n the qual- 
ity of the reddening estimates in the IHarrisl (Q996) catalogue. 
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FIG. 8. — Top panel: Ratio of our computed population-synthesis V-band 
mass-to-light ratios for LMC, SMC, and Fornax clusters (from Column 6 of 
Table 151. to those calculated by Mackev & Gilmore 1 2003a be). Our values 
are always lower because MG03 adopt a significantly steeper stellar IMF and 
do not account for gradual cluster mass loss from stellar-wind and supernova 
debris. Bottom panel: Difference of our computed population-synthesis red- 
dening minus published reddenings for individual clusters in the LMC, SMC, 
Fornax, and the Milky Way. Reddenings of Fornax and Milky Way GCs 
are taken from the compilations of Mackev & Gilmore 1 2003c ) and Harris 
i 1996) and are largely based on CMD studies of individual clusters. Redden- 
ings of LMC cluster s are co mpared to an average (E(B— V)) = 0. 1 assumed by 
Mackev & Gilmore 1 2003a), except for R136 = 30 Doradus. SMC reddenings 
are com pared to an average (E(B—V)) = 0.05 assumed by Mackev & Gilmore 
1 2003b). 



The comparisons for the old Fornax glob ulars (with true mea- 
sureme nts culled from the literature by iMackev & G ilmore 
2003c) are likewise very favorable. In the LMC and SMC 
samples, the "published" E(B-V) used for comparison in 
Fig. are generally the rou gh average values adopted by 
Mackev & Gilmore (2003a b): a constant E(B-V) = 0.10 mag 
in the LMC, and constant E(B-V) = 0.05 mag in the SMC. 
It is encouraging that our more detailed modeling returns val- 
ues clustered well around these reasonable averages. The one 
LMC cluster to which MG03 assign a non-average reddening 
is Rl 36=30 Doradus, for which they quote E(B-V) = 0.38 
mag from the CMD analysis of Hunt er et alJ {l995). Our 
model value for 30 Dor is 0.407 mag. 

In what follows, then, we make use of the population- 
synthesis E(B-V) values that we have derived for each LMC 
and SMC cluster to correct their observed surface bright- 
nesses and magnitudes for extinction. There are two LMC 
clusters (SL-663 and SL-855) for which we were unable to 
find {B - V) colors in the literature to estimate their redden- 
ings, so we simply assign to them the average (0.096) of the 
other LMC model E(B-V) values. This also happened for 
the SMC cluster NGC 361, to which we therefore assign the 
average (0.069) of the model results for the other SMC clus- 
ters. For the Fornax and Milky Way GC samples, we use the 
me asured E(B - V ) tabulated by Mackev & Gilmore (2003c) 
and lHarrisI i ll 9961) (see these papers for references to the orig- 



inal determinations of the reddenings, most of which come 
from direct study of the cluster CMDs). For every cluster in 
all four galaxies we apply our population-synthesis Ty ratios 
whenever we need to convert between luminosity and mass. 

4. DYNAMICAL MODELS AND FITS 
4.1. Dynamical Models 

We fit three type s of model to each cluster: 

First is the usual King ( 1966) single-mass, isotropic, mod- 
ified isothermal sphere, which is defined by the stellar distri- 
bution function 



f(E) oc 



exp[-£/cr^ 




]-l 



E <0 
E >0 



(1) 



where E is the stellar energy. Under certain restrictive 
conditions, this formula roughly approximates a st eady- 
state solution of the Fokker-Planck equation (e.g., iKing l 
11965b . It has, of course, already been fit to all of the 
Galactic GCs in the surfac e-brightness profile database of 
Trage r. King. & Diorgovski ( 1995), and more. These fits are 
the basis of the standard catalogues dDi orgovskil 1 993llHarrisl 
1996) of globular cluster structural and dynamical properties. 
King ( 1966) models have also been fit to a number of old 
globular clusters in other galaxies, including the LMC, and a 
few younger Magellanic Cloud clusters; but a systematic and 
uniform comparison of these standard models against all of 
the high-quality profile data that we now have for the LMC, 
SMC, and Fornax clusters has not yet been made. 

Second is a power law density profile with a core (also var- 
iously known as Moffat profiles, modified Hubble laws, or 
"Elson-Fall-Freeman" [EFF] profiles). With I(R) the lumi- 
nosity surface density of a cluster, such that fi = constant - 
2.5 log [I(R) /L Q pc~ 2 ], these models are defined by 



I(R) = 



(7~3)L t[ 
2irr 2 



l+(R/r ) 



-(7-l)/2 



(2) 



corresponding to a three-dimensional luminosity density pro- 
file j(r) oc [1 + (r / roj 1 ]' 1 1 1 . Since the density is non-zero 
even as r — ► oo, 7 > 3 is required if the integrated lumi- 
nosity is to be finite. Equation Q is strictly an ad hoc 
fitting function with no underlying basis in theory, but it 
is frequently fit to massive young clusters in the Magel- 
la nic Clouds and other galaxies, following the seminal work 
of lElson. Fall. & Freeman! ( 11987b showing that IKingl ( 119661) 
models cannot always account for the spatially extended ha- 
los of such objects. MG03 accordingly fit this model to all 
of the clusters in their sample — including the old globulars 
in the Clouds a nd the Fornax dwarf and e ven a handful of 
Milky Way GCs (M ackev & Gilmorel2003d) . However, equa- 
tion d2} has not been fit to any modern surface-brightness data 
for large numbers of Galactic globulars such as we have at 
hand. 

The physical parameters extracted for young massive clus- 
ters fr om fits of equation © by Elson, Fall, & Freeman 
( 119871) . MG03, and others (e.g.. lLarsenll2004l) are generally 
confined to the central intensity Iq, the scale radius (and 
perhaps an associated core radius R c and/or effective radius 
Rh), the power-law exponent 7, and the total luminosity L tot . 
This small set of quantities fa lls well sh ort of the wide array 
that has been calculated for King ( 1966) models, including 
relaxation times, central escape velocities, global binding en- 
ergies, and more. A full comparison between the systemic 
properties of globular and other massive star clusters would 
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seem to require comparable levels of detail in the modeling of 
all objects, no matter what specific model is actually applied; 
and a clear understanding of how model-dependent are the 
values of the derived physical parameters is critical. These 
issues can only be tackled by completely fitting both of the 
model types just described to all of the clusters in all four 
galaxies of our sample. 

Certainly a good reason why power-law models have not 
been fit systematically to Galactic globular cl usters is that 
they lack the tidal cut-off radius so important to lKind i ll 9661) 
models and obviously present in real GCs. We there- 
fore expect a priori that the power-law models may provide 
poorer fits than King models to many globular cluster den- 
sity profiles. Conversely, we already know for a fact (e.g., 
Elson, Fall, & Freeman 1987) that power laws provide better 
fits than King ( 1966) models to the outer envelopes of some 
young LMC clusters. We therefore wish also to fit all of our 
clusters with a third model which is intermediate between 
these two, o ne which can afford more extended halos than 
King ( 1966) models but still goes to zero density at a finite ra- 
dius. Another reason to c onsider such an additional mod el is 
the original suggestion of Elson, Fall, & Freeman ( 1987) (see 
also the recent review of Schweizer 2004), that the power- 
law fits to some young LMC clusters represent unbound halos 
of stars around relatively recently formed systems that could 
be stripped away over Gyr timescales, leaving behind more 
standard King-type configurations. In assessing this idea, it 
is worthwhile to ask whether self-consistent, non-King mod- 
els with extended but bound stellar halos are capable at all of 
describing these data. 

The third model that we have fit is again based on a spec- 
ified stellar di stribution function, motivated by the work of 
iWilsonl (119751) on modeling elliptical galaxies: 
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which is a different ty pe of single-m ass and isotropic modi- 
fied isothermal sphere. Wilson] i 19751) included a multiplica- 
tive term in the distribution function depending on the angular 
momentum in order to create axisymmetric model galax- 
ies. We have dropped this term from f(E) to make spheri- 
cal and isotropic cluster models, but we still refe r to equation 
as Wilson's model. The connection with the lKingl i 19661) 
model in equation Q is clear: the extra +E/oq in the first 
line of equation (|3} is simply taking away the linear term in 
the Taylor series expansion of the fundamental exp(-E /eg) 
near the zero-energy (tidal) boundary of the cluster. Although 
patently an ad hoc thing to do, the net effect of this more grad- 
ual lowering of the isothermal sphere is to produce clusters 
which are spatially more extended than King ( 1966) models, 
but still finite. 

It will be noted that there is a slight asymmetry between 
equations {1} and (|3} and equation (0: the former explicitly 
involve a velocity scale parameter <7o, while the latter incor- 
porates an explicit radial scale ro. In the formulation of his 
model. lKind l fl96 6) defined a radial scale associated with ctq: 
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where pQ is the central mass density of the model. We 
adopt the same definition for our single-mass, isotropic Wil- 
son models; and we use it also in our construction of power- 
law models to define a velocity scale ctq in terms of ro from 
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FIG. 9. — Comparison of the projected density/surface-brightness profiles 
of a single-mass King 1 1966) model cluster (dashed curves) defined by equa- 
tion Q; a spherical and isotropic Wilson 1 1975) model (solid curves) defined 
by equation J5J; and a cored power-law model (dotted curves) defined by 
equation |2|- The three examples shown all have the same dimensionless total 
luminosity Ltot/hRj r Top panel is a log-log plot of the models highlighting 
their core sturctures; bottom panel is a log-linear representation emphasizing 
their halos. Structural differences between the models are most significant 
beyond a few projected half-light (effective) radii. 



equation (|2). It is important to recognize that ro and oo are not 
equivalent, in general, to typically observed quantities such 
as a core (half -power) radius or central velocity dispersion — 
although the connections between the theoretical scales and 
these or other observables are straightforward to calculate for 
any member of these model families. 

In their dimensionless form, King and Wilson models are 
characterized by the profiles of p/po and a/ao as functions 
of r/ro. Both profiles are fully specified by the value of the 
dimensionless central potential, Wo = -0(0) ja\ > 0. In prin- 
ciple Wo can take on any real value between and oo, with the 
latter limit corresponding to a regular isothermal sphere of in- 
finite extent. Wo bears a one-to-one relationship with the more 
intuitive concentration parameter, c = log(r,/ro), where r t is 
the tidal radius of the model cluster [p(r,) = 0]. For more de- 
tails of the relations between these model parame ters, and of 
the construction of King models in general, see lKiiigl(IT963) 
or lBinn ev & Tremaine ( 1987). We have esse ntially followed 
the prescription of Binne v & Tremairiel dl987l) — with the ob- 
vious substitution of the distribution function f(E) in equation 
Q for that in equation Q — to compute Wilson models with 
arbitrary Wo or c. The result is normalized three-dimensional 
density and velocity-dispersion profiles, which are then pro- 
jected onto the plane of th e sky (using standard integ rals that 
can also be found, e.g., in lBinnev & Tremaind[T9 87) for fit- 
ting to data. 

Fixing Wo or c at some value essentially defines the over- 
all shape of the internal density profile (and the velocity- 
dispersion profile) of a King or Wilson model. The anal- 



Massive Star Clusters 



19 



ogous shape parameter for the cored power-law models of 
equation is the index 7 > 3, the exponent of the power 
law hypothesized to describe the asymptotic behavior of a 
cluster's density distribution. Clearly in this case 7 does not 
correspond to any measure of a tidal radius — r, is always in- 
finite in these models — but it still has a one-to-one connec- 
tion with the (finite) central potential of a cluster. Starting 
from equation (|2ji for the observable surface-density profile 
in these models, which is fitted directly to real cluster data, 
it is a simple matter to compute the deprojected volume den- 
sity j(r) analytically and then s olve the spherical Jeans equa- 
tion JBinnev & Tremainell 19871) . assuming unit mass-to-light 
ratio and velocity isotropy, to obtain a normalized velocity- 
dispersion profile for any 7. The complete structural and dy- 
namical details of power-law clusters are then known in full, 
and a ll of the derived physical parameters that we present in 
34.21 can be evaluated equally well within any of the three 
models that we fit. 

It is instructive first to consider how these different types of 
models compare with one another. Thus, in Figure[9]we show 
normalized surface-density profiles of one example each of a 
King, Wilson, and power-la w model cluster. We first calcu- 
lated a c = 1 .5 iKingl J1966I) model — representative of an av- 
erage Galactic globular cluster — and found its dimensionless 
total luminosity, L tot /IoRj v in terms of an arbitrary central sur- 
face density Iq and projected half-light (effective) radius R/,. 
The surface-density profile I(R), normalized by Iq and scaled 
in projected radius by R/, is shown as the dashed curves in 
both panels of Fig. [9] The top panel here is a log-log plot of 
the density vs. radius; the bottom is a log-linear plot, a repre- 
sentation which emphasizes the outer halo structure over the 
inner core regions. 

We then soug htthe lWi1sonlllT975|) model with the same di- 
mensionless L tot /IoRl as the c = 1 .5 King ( 1966) cluster. This 
turns out to be given by a Wilson c = 2.08, and the projection 
of this model onto the plane of the sky is shown as the solid 
curves in Fig.|9] Similarly, the 7 = 3.765 power-law model 
shown as the dotted curves has the same L tot /IoRj 7 again. 
Thus, Fig.|9]illustrates the relative spatial extent predicted by 
these three models for a cluster with fixed total luminosity, 
central surface brightness, and projected half-light radius (all 
quantities that tend to be observationally well-determined in 
real clusters). Clearly apparent is the intermediacy of the Wil- 
son model between the more sharply truncated King model 
and the infinite power law. Also evident is that the differences 
between these models are generally largest beyond a few ef- 
fective radii, in the outer halos of clusters. 

Fi gure ITol extends this comparison to general Kin gl (119661) 
cluster concentrations. [A more comprehensive discussion of 
the defining features of King models v s. Wilson models in 
particular can be found in Hunter ( 1977).] Here, we have cal- 
culated the dimensionless L tot /IoRj r for each of a large number 
of King models with 0.3 < c < 4, and then found the unique 
Wilson c and power-law 7 that give the same dimensionless 
luminosity. Although we have done this formally without any 
restrictions on the parameter c, note that the majority of King- 
model fits to real GCs return c < 2, and indeed the model itself 
is unstable to the gravothermal catastrophe at concentrations 
higher than this. 

It is noteworthy from this figure that, first, tow- 
concentration King or Wilson models, which are character- 
ized by a sharp decline in density beyond a dominant, nearly 
constant-density core, find their analogue in high-j power 
laws. Second, while the shape parameters c or 7 are certainly 
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FIG. 10. — Relation between King-model concentration c = log (r f /ro), 
Wilson-model c, and power-law exponent 7 when a cluster is required to have 
the same total luminosity and effective radius and central surface-brightness 
in all three models. (Each c or 7 corresponds to a unique value of L\at/IoRi.} 
King and Wilson models with c < exist in principle (c — ► —00 as the dimen- 
sionless central potential Wo — » 0) but are not shown here, as real clusters of 
such low concentration are rare. The same is true of power-law models with 
7 > 10. 



useful indicators of a cluster's global structure, their numeri- 
cal values are highly model-dependent. This makes it of in- 
terest to consider a more generally applicable "concentration 
index" able to represent the spatial extent or potential depth 
of any cluster in a more model-independent way (so as to al- 
low, e.g., for a combined analysis of clusters which may not 
all be described well b y the same type of model). We return 
to this point below in 34.21 where we now present the fits of 
all models to our cluster sample. 

4.2. Surface-Brightness Fits and Cluster Physical 
Parameters 

Our fitting procedure involves computing in full large num- 
bers of King, Wilson, and power-law structural/dynamical 
models, spanning a wide range of fixed values of the appro- 
priate shape parameter Wo or 7. Separately for each fam- 
ily in turn, we then fit every model on the appropriate Wo- 
ol 7-grid to a cluster's observed surface-brightness profile, 
Hv = fiy.Q-2.5 log [I(R/ro)/Io], finding the radial scale ro and 
central SB /iy.o which minimize x 2 f° r every given value of 
Wo or 7. The (Wo> r O,^v,o) or (7j r O ) A t v,o) combination that 
yields the global minimum x^in over tne SX^ use d defines the 
best-fit model of that type. Estimates of the one-sigma uncer- 
tainties on these basic fit parameters, and those on all associ- 
ated derived quantities, follow from their extreme values over 
the subgrid of fits with x 2 < Xmm + 1 ■ F° r tne most P art ' we 
calculate and minimize \ 2 as the weighted sum of squared dif- 
ferences between model and observed intensities / (in units of 
V-band L Q pc~ 2 ), rather than logarithmic surface-brightness 
units; although in a few cases more stable fits resulted from 
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defining x 2 as the weighted sum of squared surface-brightness 
deviations. The vast majority of our fits are error-weighted, 
with the uncertainties on individual datapoints either taken 
from the original sources of HST and ground-based starcounts 
in the LMC, SMC, and Fornax clusters, or es tima ted by us 
from the Trager et al. (1995) catalogue data ( 32.21 and Table 

B 

4.2.1. Fits to lj Centauri 

As a detailed example, Figure ^2 shows our best King, 
Wilson, and power-law fits to the well-studied Galac- 
tic globular cluster lj Cen tauri = NGC 5139 (see also 
McLaugh lin & Mevlanl2 003). The main parameters and min- 
imum x 2 values for the three fits are listed in Table [9] These 
are also given, along with many other derived structural and 
dynamical parameters, in our results below for the full cluster 
sample (Tables[Tolthrough ll4t . 

The upper panels of Figure ^2 display the cluster's model 
and observed surface-brightness profiles, while the lower 
panels show the predicted and observed internal velocity- 
dispersion profiles, which we consider in detail for this one 
object onl y. The V-band surface-brightnes s datapoints are 
those from Tra ger. King . & Diorgo vskil ill 9951) with our esti- 
mated errorbars attached. The radial scale is given in units of 
arcseconds along the lower horizontal axes and parsecs (for an 
assumed distance of 5.3 kpc taken from Harris 1996) along the 
upper axes. The left-hand upper panel is a log-log representa- 
tion of density vs. radius, which focuses on the core structure 
of the cluster; the right-hand panel is in log-linear format to 
highlight its halo structure. 

The main point to be taken from the top of F ig. II His that 
even among Milky Way globular clusters, where King ( 1966) 
models are generally taken as the best physical description, 
alternate models can fit even the most comprehensive data 
at least as well and sometimes, as in the case of the Wilson 
model here, better. To be sure, the single-mass and isotropic 
King models that we fit have long been known to be inade- 
quate for many old GCs, a fact that is usually attributed by 
hypothesis to the influences of energy equipartition and mass 
segregation and/or velocity anisotropy in well-relaxed mul- 
timass stellar populations. Multimass and anisotropic mod- 
els still based on the King distributi on function (eq . ffl ) 
have therefore been developed ( Da Costa & Freemar] ll976t 
iGunn & Griffinin"979ft and fit to a go od number of Galac tic 
globulars — including lj Centauri (e.g., Mevlan et al. 1995) — 
— with much better success. The case of lj Cen is instructive, 
however, in that it is an unrelax ed cluster showin g no evidence 
of advanced mass segregation ( Anderson 199 3), nor of veloc - 
ity anisotropy jMerritt et alJll997t Ivan Leeuwen et al1l2000t) . 
Compensating for the shortcomings of the simple lKma (1966) 
model in Fig.^Jis therefore better done with entirely differ- 
ent single-mass and isotropic models — even if they are ad hoc 
to some degree — than with multimass and anisotropic varia- 
tions on the King distribution function. This is an additional 
justification for our fitting of Wilson and power-law models 
to a full suite of Galactic globulars. lj Cen is the most mas- 
sive GC in the Milky Way, and an uncharacteristically diffuse 
one. It might therefore have been thought, justifiably, to be 
a rather special case; but as we shall see, our single-mass, 
isotropic Wilson models in particular do provide fits of qual- 
ity comparable to or better than King models for a majority of 
GC surface-brightness profiles. We also note that the cluster 
tidal radius implied by the Wilson-model fit in lj Cen is only 
50% larger than that of the King model (see Tablell llbelow). 



meaning that the generally greater extent of Wilson's model 
relative to King's need not imply severe inconsistencies with 
the expected sizes of GCs for a given Galactic tidal field. 

The bottom panels of Fig. II II show the obs erved velocity 
dispersion as a fun ction of radius in lj Cen llMevlan et alJ 
1995; Seitzer 1983), again with a logarithmic radial scale on 
the left and a linear radial scale on the right. The model 
curves now are the profiles predicted (after projection along 
the line of sight) by solving the spherical Jeans equation with 
the de-projections of the best-fit density profiles in the upper 
panels. The model predictions are inherently dimensionless 
in form, yielding a p /ctq as a function of projected radius R 
where ao is the velocity scale definedby ctq = (AnGTy jorfy/9 
(cf. eq. |4l). Knowing the radial scale rg and the central three- 
dimensional luminosity density jo from fitting to the cluster 
surface brightness, the theoretical velocity-dispersion curves 
are scaled to match the data essentially by fixing a (spatially 
constant) mass-to-light ratio Ty. The value of Ty that gives 
the best agreement with the velocity dispersions in lj Cen is 
listed for each model type in the lower right-hand panel of 
Fig. The main conclusion is that our different surface- 
brightness fits are, in general, associated with self-consistent 
internal dynamics that are closely similar both in relative 
terms (the shapes of the model velocity-dispersion profiles) 
and in an absolute sense (the implied value of Ty). 

4.2.2. Fits to All Clusters 

Figure^] which can be found at the end of the paper, dis- 
plays the best-fit King, Wilson, and power-law models for the 
V-band surface-brightness profiles of our 53 LMC clusters, 10 
SMC clusters, the 5 globulars in Fornax, and 84 Galactic GCs 
besides lj Centa uri (recall that we do not present fi ts to the en- 
tire database of Tr ag er. King. & Diorgovskil 19951 for reasons 
discussed in 32. 2\ . Each cluster is presented in two panels 
formatted as in the top of Fig. ^2 one with both axes on log- 
arithmic scales and one in log-linear form. Along the top of 
every panel, projected radius is given in pc, obtained by as- 
suming a distance of 50. 1 kpc to all clusters in the LMC, 60.0 
kpc to all clusters in the SMC, and 137 kpc to Fornax; helio- 
centric distances to indiv idual Milky Way GCs are taken from 
the catalogue of lHarrisI (119961) . In all cases, dashed curves 
are the King-model fits; solid curves are Wilson models; and 
dotted lines are power laws with cores. 

In the LMC, SMC, and Fornax cluster panels, the open 
circles correspond to su rface-brightness data taken from 
Mackev & Gilmore (2003a be), after applying the zeropoint 
corrections discussed in 32.11 (see Table |4} and correcting 
for the V-band extinction A v = 3 . 1 [(B - V) - (B - V) ] inferred 
from our population-synthesis modeling in 33.2l (see Table[8j. 
There are generally two such points at each radius, one com- 
ing from the primary V-band counts of MG03 and the other 
from their secondary B or / profiles shifted by the color terms 
in Table [2] The open squares in the plots for LMC and SMC 
clusters denote the ground-based starcount data that we have 
collected and matched onto the re-zeropointed MG03 profiles. 
Asterisks refer to datapoints that we have not included in cal- 
culating and minimizing x 2 to identify the best model fits (i.e., 
those with "fit flags" of in Table[5j. 

In the Mil ky Way GC panels, the open circle s are the data- 
points from lTrager. King. & Diorgovskil(ll995l) . taken exactly 
as published except for Palomar 10 and Terzan 7, which we 
have calibrated as described in 32.21 The errorbars on these 
points are our own estimates based on Table |6] The asterisks 
in these cases are surface-brightness measurements also pub- 
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TABLE 9 

Model Fits to w Centauri = NGC 5139 



Model 

(1) 


Reduced v 2 a 
(2) 


(3) 


c = \og(r,/r Q ) 
(4) 


Mv,o b 
(5) 


r [sec] 
(6) 


D [pc] c 
(V) 


King 

Wilson 

Power-law 


1.73 
0.47 
1.47 


W = 6.20+°, 2 ° 
Wo = 4.70t|» 
^-5 55+0.10 


1 31+0.05 

1 34+8:83 

1,J -0.03 
OO 


16.44+° 06 
16 52+°' : °^ 
16 65+° 05 


141.20+"° 
196.81!"? 
327.58+H^ 


3.63+0" 
5 ofi+ttii 

• ,UD -0.18 
8 42+ ' 29 



a There are 54 datapoints in the fitted profile, and we fit by finding the fiy.o antl r o which minimizes \ 2 over a large grid 
of fixed Wq or 7 values. Thus, the reduced Xmin P er d e g ree of freedom is just Xmi n /52. 

^Corrected for extinction given in Harris 1 199d) . 

c Assumes a heliocentric distance of D = 5.3 kpc IHarriil996l) . 
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FIG. 1 1 . — Detailed fitting of Kina 1 1966), Wilson 1 1975), and power-law models to the Galactic globular cluster ui Centauri. Top panels show the surface- 
brightness profile as a function of logarithmic projected radius on the left, linear projected radius on the right. Bottom panels show the observed internal 
line-of-sight velocity-dispersion profile from Mevlan et al. 1 1995, open squares) and Seitzei 1 1983, filled triangles), and the model profiles after normalization by 
the V-band mass-to-light ratios indicated. 
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lished by Trager at al. but given relative weightings w, < 0.15 
by them; we did not include these points when fitting our 
models. 

Table 1 101 (also at the end of the paper) gives a few impor- 
tant numbers for each cluster and lists the main parameters 
of every fit. The first column here is the cluster name with 
a short prefix identifying its parent galaxy. Column (2) is 
the zeropoint shift that we have applied to published surface- 
brightness magnitudes (copied directly from Table @] for the 
MG03 cluster sample; for all but Palomar 10 and Terzan 7 
in the Milky Way GC sample). Column (3 ) is the estimated 
V-band extinction, derived as discussed in 33.21 Next are the 
assumed distance to the cluster, the number of datapoints ulti- 
mately included in fitting any model to the cluster, and a code 
for the weighting scheme used to minimize \ 2 - 

Subsequent columns in Table 1 101 cover three lines for each 
cluster, one line for each type of model fit. Column (7) iden- 
tifies the model; Column (8) gives the minimum unreduced 
X 2 obtained for that class of model; Column (9) gives the 
appropriate "shape" parameter Wo or 7 at which x 2 is min- 
imized; Column (10) gives the concentration c = log(r f /ro) of 
the best fit (this is related uniquely to Wo for King and Wil- 
son models but is always 00 for power-law models); Column 

(11) gives the best-fit central surface brightness after correc- 
tion for both extinction and any zeropoint change; Column 

(12) gives the best-fit model scale radius, ro, in arcseconds; 
and Column (13) gives the value of ro in parsecs. The un- 
certainties in Wo, c, 7, and ro reflect their variations among 
model fits with x 2 < Xmm + 1 • ^ e enorbars on jiyfi derive 
from these formal fitting uncertainties combined in quadra- 
ture with the uncertainties in the zeropoint offsets Apy and 
extinctions Ay (the x 2 of any fit is calculated using the pub- 
lished uncertainties in the individual SB datapoints before any 
systematic corrections are applied). 

Tablell ll( at the end of the paper) contains a number of other 
structural cluster properties derived from the basic fit param- 
eters: 

• log r, = c + log ro is the model tidal radius. 

• log R c refers to the projected core radius of the model 
fitting a cluster. It is defined by I(R C ) = h/2, or 
Hv(Rc) — Mv.o + 0.753, and is not the same in general 
as the radial scale ro in Table 1 101 For King and Wilson 
models, ro is simply a convenience defined by dimen- 
sional analysis of Poisson's equation, although it does 
bear a uniqe relationship to the observable R c through 
well-defined functions of Wo or c. For power-law mod- 
els, the connection between ro and R c stems from the 
defining equation (|2j above. 

• log R i t refers to the half-light, or effective, radius of a 
model: that radius containing half the total luminosity 
in projection. It is related to ro by one-to-one functions 
of Wo or 7. 

• log(Ri t /R c ) is a measure of cluster concentration that is 
relatively more model-independent than Wo or c or 7, 
in the sense that it is generically well defined and its 
physical meaning is always the same. We consider it a 
more suitable quantity to use when intercomparing the 
overall properties of clusters which may not all be fit 
by the same kind of model (cf. our earlier discussion 
around Fig.lTol. 



• log Iq = 0.4 (26 .422- /1 v/o) is the logarithm of the best- 
fit central luminosity surface density in the V band, 
in units of L Q pc~ 2 . The surface-brightness zeropoint 
of 26.422 correspo nds to a s olar absolute magnitude 
M v . e = +4.85 fe. g ..lT. a nd1999h. 

• log y'o is the logarithmic central luminosity volume den- 
sity in the V band, in units of L Q pc~ 3 . It is given by 
jo = Jhjro, where J' is a smooth, model-dependent 
function of Wo or 7, which we have calculated in detail 
for King, Wilson, and power-law models. 

• log L tot is the logarithm of the total integrated model 
luminosity in the V band. It is related to the product 
Io r o by model-dependent functions of Wo or 7. 

• y tot = 4.85-2.5 log(Lt ot /L ) + 5 log(D/10pc) is the to- 
tal, extinction-corrected apparent magnitude of a model 
cluster. 

• log//, = log(L tot /2nRj r ) is the V-band luminosity sur- 
face density averaged over the half-light/effective ra- 
dius, in units of L pc" 2 . 

• {nv)h = 26.42-2.5 log(//,/Lopc~ 2 ) is the average sur- 
face brightness inside the half-light radius, in V-band 
mag arcsec" 2 . 

Again, the uncertainties on all of these derived parameters 
have been estimated (separately for each given model family) 
by calculating them in every model which yields x 2 within 
1 of the global minimum for a cluster, and then taking the 
differences between the extreme and best-fit values of the pa- 
rameters. 

Table [21 (see the end of the paper) next lists a number of 
cluster properties derived from the structural parameters al- 
ready given plus a mass-to-light ratio. The first column of 
this table contains the cluster name, as usual. Column (2) 
shows the mass-to-light ratio, in solar units, that we have 
adopted for each object from the analysis in 5j3] — that is, on 
the basis of population-synthesis modeling given individual 
ag es and metallicities for th e clusters, using the model code 
of Bruzual & Chariot (2003) and assuming the disk-star IMF 
of Chabriei (2003). The values of Y pop in Tablell2lhave been 
copied directly from Column (6) of Table [8] The remaining 
entries in Tablell2lare. for each type of model fit to each clus- 
ter: 

• log M tot = log Y pop + log L tot , the integrated model mass 
in solar units. 

• log Eb, the integrated binding energy in ergs, defined 
through Eb = -(1/2) f r ' 4nr 2 p(j)dr. Here the minus 
sign makes Eb positive for gravitationally bound ob- 
jects, and (f>(r) is the potential generated (through Pois- 
son's equation) by the model mass-density distribution 
p(r). Ei, can be written in terms of the fitted central lu- 
minosity density jo, scale radius ro, a model-dependent 
function of Wo or 7, and Y pop . A more detailed out- 
line of this procedu re for King models may be found in 
McLaughlin (2000), which we have followed closely to 
evaluate Eb for our Wilson and power-law fits as well. 

• log So = log Y pop + log /o, the central mass surface den- 
sity of the model fit in M Q pc~ 2 . 
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• logpo = log T pop + log jo, the central mass volume den- 
sity in Mq pc~ 3 . 

• log £/, = log Ty° p + log //,, the model mass density aver- 
aged over the half-light radius R/, (which is equal to the 
half-mass radius under our assumption of single-mass 
stellar populations, i.e., spatially constant mass-to-light 
ratios). 

• <t Pi o, the predicted line-of-sight velocity dispersion at 
the cluster center, in km s" 1 . As was already suggested 
above, the solution of Poisson's and Jeans' equations 
for any model yields a dimensionless ct p ,o/cto, and with 
(To given by the fitted ro and p through equation @, the 
predicted observable dispersion follows immediately. 

• v esc ,o, tne predicted central "escape" velocity in km s" 1 . 
A star moving out from the center of a cluster with 
speed v eS c ,o will j ust come to rest at infinity. In general, 
then, v^ sc /cr^ = 2 [\V + GM tol /r,al] . Note that the sec- 
ond term on the right-hand side of this definition van- 
ishes for power-law models, in which r, — > oo. In these 
models a (finite) dimensionless Wo is associated with 
every value of 7 > 3 by solving Poisson's equation with 
0(oo) = O. 

• log f r h, the two-body relaxation time at the model 
projected half-mass radius. This is estimated 

as f r h/yr = [2.06 x 10 6 /ln (0.4M, nt /OTj] rn^M^R^ 2 
(Binn ev & Tremainell 19871 eq. 8-72), if (the aver- 
age stellar mass in a cluster) and M lot are both in solar 
units and /?/, is in pc. We have evaluated this timescale 
assuming an average m* = O.5M in all clusters. 

• log /o = log [po / (27rtr^) 3 / 2 ] , a measure of the 
model's central phase-space density in units of 
M Q pc" 3 (kms -1 ) -3 . In this expression, a c refers to the 
central one-dimensional velocity dispersion without 
projection along the line of sight. The ratio a c /ao is 
obtained in general from the solution of the Poisson 
and Jeans equations for given Wo or 7, and the fitted 
do again is known from equation ©. With the central 
relaxation time t rc of a cluster def ined as in equation 
(8-71) of Binnev & Tremaine ( 1987), taking an average 
stellar mass of m* = 0.5M Q and a typical Coulomb 
logarithm In A w 12 leads to the approximate relation 
log(f rc /yr)~8.28-log/ . 

The uncertainties in these derived dynamical quantities are 
estimated from their variations around the minimum of \ 2 on 
the model grids we fit, as above, combined in quadrature with 
the population-synthesis model uncertainties in T^ p , which 
in turn reflect the estimated uncertainties in published cluster 
ages and metallicities. 

4.3. Fit Comparisons 

4.3.1. Our Fits vs. Published Catalogues 

An imp ortan t check on the array of cluster properties pre- 
sented in 34.2l is provided by comparing our basic fit parame- 
ters against those in published catalogues for the Milky Way 
and LMC/SMC/Fornax clusters. We confine such compar- 
isons to one betwee n our power- law fit p arameters and those 
reported bv lMackev & Gilmord fe003albld MG03) for LMC, 
SMC, and Fornax clusters; and one between our King-model 



fit s to the Galactic GCs and the parameters in the catalogue 
of lHarrisI (|1996) (its elf ultimately based in large part on the 
profile database of iTrager. King. & Diorgovskil [l995). The 
generally good results give us confidence that there are no 
problematic biases or any other procedural issues with our 
modeling or fitting techniques; thus, we have not searched 
the literature to compare with all other model fits that might 
have been performed on any of these clusters. 

The left-hand panels of Fig.ll3lshow our power-law expo- 
nents 7, scale radii ro, and central surface brightnesses py,o 
against those published by MG03 for their cluster sample. 
The overall agreement between our numbers and theirs is 
apparent — although note that the "Mackey & Gilmore" sur- 
face brightnesses we compare to in the bottom left panel have 
already been corrected for our extinctions and zeropoint off- 
sets in Table 1 101 There is somewhat more scatter in the plots 
of 7 and ro values than in the nv,o g ra ph> for the simple rea- 
son that the latter parameters are more sensitive to the ground- 
based starcounts which we have added to about two-thirds of 
the MG03 cluster sample. 

It is worth noting a frequent tendency for our power-law fits 
to return somewhat steeper 7, and correspondingly larger ro, 
than those of MG03. This can also be traced to our inclusion 
of ground-based data, as the specific example of the LMC 
cluster NGC 2121 illustrates well. MG03 quote a value of 
7 = 3.25 for this object, whereas we have obtained 7 = 6.00; it 
is clearly visible as an "outlier" in the upper left-hand panel of 
Fig-El Referring back to the plotted fits in Fig.[l2] however, 
we see that the SB profile from Mackev & Gilmore (2003a) 
extends only to R ~ 70" = 17 pc in this case — just barely out- 
side the constant-density core region of NGC 2121 — while 
our additional groundbased data reach to R ~ 100" ~ 24 pc 
and are critical to accurately constraining any model fit. For- 
tunately this example is the most extreme in our sample, but 
the point is made that when a model extrapolates significantly 
beyond the limit of the fitted data, care must be taken in using 
the results. This is of most concern for power-law models, 
which are innately the most spatially extensive of the ones we 
fit. 

The right-hand panels of Fig. ^] show our fitted King- 
model concentrations, scale radii, and central surface bright- 
nesses for 85 Galactic globular clusters, against those param- 
eters taken from the Harris catalogue. For the most part, the 
agreement is again quite good; but there are a number of clus- 
ters for which we claim s ignificantly different c values from 
those indicated by Harris ( 1996), and these differences prop- 
agate into the ro and fiy,o plots (given the same data, a lower 
fitted c corresponds in general to a measurably larger ro and 
somewhat fainter fivo)- 

For the two GCs NGC 6101 and NGC 6496, we find King 
concentrations c«0 while Harris gives c ~ 0.7-0.8. Look- 
ing at these objects in Fig. ^] they are unquestionably low- 
concentration clusters, with steep declines in surface bright- 
ness beyond relatively large cores. The exact details of any fit 
are heavily influenced by the outermost one or two datapoints 
in each case (and by irregular structure at the center of NGC 
6496), and our errorbars on c (and ro and py,o) reflect this. 

Some what similarly , we find c = 1 .5 ± 0. 12 for NGC 6528, 
whereas Harris ( 1996) sta tes c = 2 . 3. Th ere are many points 
from [Trager, King, & DiorgovskH i 19951) for this cluster to 
which we assigned zero weight when doing our fits. The 
model parameters are apparently quite dependent on which 
of the data are taken into account for this relatively poorly 
defined brightness profile. 
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FIG. 13. — Comparison of our power-law fit parameters for LMC, SMC, 
and Fornax clusters (left-hand panels) against those from Mackev & Gilmore 
i 2003a b c), and of our King-model fit parameters for Milky Way globular 
clusters (right-hand panels) against those catalogued by Harris 1 1996). Bro- 
ken lines in all panels indicate equality. The comparison of our central sur- 
face brightnesses with the "Mackey & Gilmore" values in the lower left-hand 
panel uses their dat a aft er correction for both the zeropoint offsets Afiy m 
Column (2) of T able [Tol (also Tablel4l and the V-band extinctions in Column 
(3) of Table ITol The Milky Way GC central surface brig htne sses plotted in 
the bottom righ t panel are co rrected for the Ay in Table ITol which are the 
same as those in lHamHT9^6T) . 



The three points falling highest above the line of equality 
in the top right-hand panel of Fig.ll3lcorrespond to Palomar 
10 (Harris c = 0.58; our c = 1.59), Palomar 1 (Harris c = 1.6; 
our c = 2.57), and Palomar 12 (Harris c = 1.94; our c = 2.98). 
Inspection of these in Fig.^lshows no obvious problem with 
any of our models f or Pal 10, al though clearly th e range of fit- 
ted datapoints from Trager, King, & Diorgovski ( 1995) is less 
than ideal; in fa ct, Harris' tabulated parameters are ba sed on 
observations by Kaisler, Harris, & McLaughlin ( 1997) which 
supersede those of Trager et al. The data for Pal 1 do not show 
any clear evidence for an isothermal core in the first place, 
allowing this cluster to be fit by a high-concentration King 
model with a small predicted scale radius and bright central 
surface brightness that are in fact unobserved. Pal 12 shows 
what might be described as a double-core structure, and it is 
therefore also fit relatively best by high-c, low-ro and bright- 
fiv.o King and Wilson models — although the fits are certainly 
not "good" in an absolute \ 2 sense (see Table II Oi l. Any cata- 
logued fit parameters (including ours) for these three Palomar 
clusters are probably best viewed as only provisional. 

4.3.2. Goodness-of-Fit for Different Models 

We next compare the \ 2 values of the different model fits to 
every cluster in our sample. We make this comparison relative 
to x 2 °f the best-fit King model in each case by computing 
A = (x 2 -Xking)/(X 2 + X K in g ) for the best Wilson and power- 
law fits. This index ranges from a minimum of A = -1 for an 
alternate model with x "C Xking' t0 a maximum of A = +1 



when xling < X 2 ■ 

The top panels of Fig.H4lshow the distribution of A values 
for the Wilson (filled circles) and power-law (open squares) 
fits to the MG03 cluster sample in the LMC, SMC, and For- 
nax. The points for the Wilson and power-law fits of any one 
cluster are connected by a solid line. On the left we plot A 
as a function of cluster dynamical age, i.e., the chronological 
age r in units of the King-model half-mass relaxation time 
f,h from Table [21 On the right, we show A as a function of 
the ratio Ru st /Rh, where /?i ast is the clustercentric radius of the 
outermost surface-brighntess datapoint observed in a cluster 
and R/, is the (King-model) half-mass radius. 

It is immediately apparent that, in every LMC/SMC/Fornax 
cl uster s t udied here, iWilsonl (11975) models fit at least as well 
as iKingl (11966) models, and very often substantially better. 
Power-law models generally fit roughly as well as Wilson 
models, sometimes slightly better and sometimes somewhat 
worse. On one level, the top panels of Fig.^Jare therefore a 
re-statement of the appreciated fact that many clusters in the 
Magellanic Clouds are more ex tended tha n classic King mod- 
els (e.g. , [Elson. Fal l. & Freemanl["l987l iMackev & Gilmorel 
l2003albl) . New here is the demonstration that, although 
they are usually acceptable fits, untruncated power-law forms 
specifically are not required to describe the density distribu- 
tions of these objects. We expect that the same is likely true 
of young massive clusters in other disk galaxies ( iLarsenl 2004; 
ISchweizeJ2004l and references therein). 

Also new is our quantification of the improvement in fit 
yielded by the more extended models as a function of clus- 
ter age. Elson, Fall, & Freeman ( 1987) originally suggested 
that power-law models fit young LMC clusters better than 
King models because of the presence of unbound stellar halos 
which are relics of the cluster formation process and simply 
have not had time to be stripped away by tides. The chrono- 
logica lly young clusters in the curren t sample (which includes 
that of Elson, Fall, & Freeman 1987) generally have r/f r h < 1 
in the upper left panel of Fig.^J Many of them have A < 
and are undoubtedly fit better by Wilson or power-law mod- 
els than by King models; but about as many have A«0 and 
are equally well fit by any model type. Moreover, there is a 
comparable number of chronologically and dynamically old 
clusters (including the Magellanic Cloud globulars and those 
in Fornax) which, at r/f r h ~ 10, also have A < and are char- 
acterized by spatially extended halos not easily reproduced by 
King models. [Working f rom completely independent data, 
Rodgers & Roberts ( 1994) have already claimed this for three 
of the five globular clusters in Fornax.] Age appears not to 
be the main factor in determining whether or not any of these 
clusters can be well described by regular King models. 

Instead, the upper right-hand panel of Fig.[l4]shows simply 
that King models provide progressively less satisfactory fits 
as surface-brightness and starcount data extend farther and 
farther into cluster halos. All of the clusters here — young 
or old, in the disk or halo populations of any of these three 
small galaxies — are better fit by Wilson or power-law models 
if the models are forced to fit beyond ~ 4-5 half-light radii in 
the clusters. (It will be recalled from Fig. ^0] above that this 
is the point where the King, Wilson, and power-law model 
structures begin to differ appreciably.) We suggest that the 
halos of massive star clusters are gener ically more e xtensive 
than the stellar distribution function of King ( 1966) allows, 
and that the development of a physically motivated model ac- 
counting for this (one less ad hoc than a Wilson or power-law 



Massive Star Clusters 



25 



1 I — i i i mil — i i 1 1 inij — i i 1 1 inij — i i i nil — i rn 1 — i — i i i 1 1 1 1 1 — i — i i i i 1 1 1 

LMC, SMC, Fornax 




0.1 1 10 100 1 10 



Vt rh R,ast/R h 

FIG. 14. — Goodness of fit of Wilson 1 1975) spheres and power-law models, relative to standard Kina 1 1966) models, for 68 clusters in the LMC, SMC, 
and Fornax (top panels) and for 85 globular clusters in the Milky Way (bottom panels). Solid points mark the relative \ 2 index, A, for the Wilson-model fits; 
open squares denote A for the power-law fits. Solid lines connect the Wilson and power-law A values for each cluster. A is shown as a function of cluster 
dynamical age in the left-hand panels, and as a function of the radial extent of the observed surface-brightness profiles in the right-hand panels. The inherently 
more extended Wilson model tends to fit clusters of any age better (A < 0) when their outer halos are better defined empirically. 



prescription) could lend substantial new insight into questions 
of cluster formation and evolution. 

The fun damental p hysical question remaining from 
Elson . Fall. & Fr eeman (1987) is whether the stars at the 
largest observed radii in the young clusters particularly are 
gravitationally bound in a model-independent sense: do these 
objects overflow the Roche lobes defined by the potentials of 
their parent galaxies? Simply fitting power laws to the clusters 
does not in fact address this issue, as equally good or better 
fits of spatially limited Wilson models can be found; but at the 
same time, it remains to be shown that the fitted r t from the 
Wilson models actually correspond to the true tidal limits im- 
posed by the galaxies. Moreover, the question now has to be 
extended to many old globulars. To properly answer it for any 
cluster requires at a minimum not only a highly precise empir- 
ical estimate of r, itself, but also a detailed understanding of 



the total mass distribution and gravitational field of the parent 
galaxy; good information on the present-day galactocentric 
position of the cluster; and knowledge of its orbital energy and 
pericenter. The interplay between these ingredients makes for 
a subtle problem, fraught with uncertainty even in the Milky 
Way (see, e.g., Ilnnanen. Harris. cC Webbink 19831) . which is 
beyond the scope of this paper. However, we do touch briefly 
on a comparison between the fitted "tidal" radii in the King 
and Wilson models for our cluster sample, in Fig.ll8lbelow. 

The bottom panels of Fig.^Jare plots of the relative \ 2 in- 
dex A against both dynamical age and spatial extent of the SB 
data for the 85 Galactic GCs that we have modeled. 4 It is clear 

4 Note the absence of Milky Way globulars with ffi ast /Rj < 1 in Fig.lMl 
contrasting with the presence o f such objects in the LMC/SMC/Fornax clus- 
ter sample. As we discussed in 82.21 such large fitted R;, in the Galactic sam- 
ple tend to be associated either with relatively poor data or with clusters iden- 
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that many globulars with well-observed halos (R\- dst /Rh > 5) 
are again relatively better fit by Wilson models than by single- 
mass, isotropic King models, regardless of their dynamical 
age r/f r h- In most of these cases power laws are worse fits 
than Wilson models, which likely reflects the inability of the 
former to describe tidal limits to a cluster. But now there 
are also about ten globulars with data extending beyond 5R/„ 
which are better fit by King models than either Wilson or 
power-law spheres. A good example is NGC 104 = 47 Tu- 
canae (A = +0.5 for the Wilson model fit), although in Fig. 1121 
this profile actually appears to prefer some kind of description 
intermediate to King ( 1966) and Wilson H 19751) . 

The bottom left panel in this plot gives the visual im- 
pression that — unlike in the LMC, SMC, and Fornax cluster 
sample — there may be some correlation between our A statis- 
tic and the dynamical age r/t^. A test using the Spearman 
rank correlation coefficient for the Wilson-model fits specif- 
ically shows that these quantities are indeed correlated, with 
a formal confidence level of > 99%. However, the correla- 
tion between A and Ri a st/Rh for the Milky Way GCs is still 
stronger and more significant than any correlation with r/f r h. 
In any case, all of this is going on over a narrow range of 
extreme age relative to the rather larger spread in the LMC, 
SMC, and Fornax sample, where the situation is much clearer. 

Some of these differences in the bottom panels of Fig. [T4l 
relative to the upper panels may simply be a reflection of 
the more heterogeneous nature of the original data compiled 
by iTrager. King. & D iorgovskil (I 19951) . It could also be in- 
dicative of the ultimate limitations of our assumptions of ve- 
locity isotropy and (probably more important) a single-mass 
stellar population in the clusters: NGC 104, for exampl e, is 
known to exhibit mass segregation (e.g., Anderson 1997) and 
can b e fit well with multimass King models (Mevlanl ll988l 
1989). On the other hand, we find here that the cluster NGC 
5272 = M 3, which was the original motivation for the de- 
velopment of multimass Kin g models (iDa Costa & Freemanl 
119761 Idunn & Griffinlll979l) . can be perfectly well fit by a 
single-mass and isotropic Wilson model; it is the filled circle 
at R\dst/Rh = 12.4 and A = -0.87 in the lower right panel of 
Fig.nl 

The totality of the results present ed here still s uggest to us 
that a fundamental alteration to the lKingl (119 66) distribution 
function is required to account for the halo structure and dy- 
namics of massive star clusters in general. That the ad hoc, 
single-mass, and isotropic l Wilson! d!975h model is not the per- 
fect solution should come as n o surpr ise; but, as our discus- 
sion of uj Centauri concluded ( 34.2. 1> . neithe r can mul timass 
and/or anisotropic variations on the standard King ( 196Q) dis- 
tribution function correctly explain the structure of all globu- 
lar clusters. 

This issue aside, we now turn to consider whether the ob- 
servable, physical properties — core and half-light radii, to- 
tal luminosities, and the like — that we have derived for our 
clusters are reasonably model-independent. In particular, we 
would like some assurance that they are robust enough to al- 
low useful characterizations of parameter interdependences 
and trends between clusters. 



4.3.3. Physical Cluster Properties in Different Models 

tified as core-collapsed by Trager, Kins. & Diorg ovsk i 1 1995). These G Cs 
are named in TablelTlbut have been left out of Fig. ll4l (and Tables [Tolll2l by 
construction. 



Figure [21 compares the Wilson- and King-model values 
for a number of parameters of the LMC, SMC, and For- 
nax (MG03) clusters in our sample (see Tables IToTll2l . The 
overriding conclusion to be drawn is that these two types 
of model fits tend, with some understandable exceptions, to 
agree within <~ 5% (<~ 0.02 dex) on the values of basic cluster 
properties. We have plotted the differences in fitted and de- 
rived properties between the two models as functions of the 
relative X 2 index A = (x^i son -xLg)/(Xwason + xLg)- 

The upper left-hand panel of Fig. [H] shows that in most 
cases the central surface brightness /i^o is very well deter- 
mined, stable at the ~ 0.03-mag level on average no matter 
which model is fit. (Note, however, that when the Wilson 
model fits better, it tends to return a slightly fainter central 
brighntess than a King-model fit.) One of the most obvious 
exceptions is the LMC cluster NGC 2005, which is the point 
enclosed in an open square; it in fact lies off the vertical scale 
at ^t/ (W)-^y,o(K) = -5.01. Although the \ 2 of the two fits 
are essentially the same, the Wilson model for this cluster is 
of much higher concentration — and therefore brighter central 
surface brightness — than the King model, due to the influ- 
ence of a single datapoint at R ~ 1" ~ 0.24 pc (see Fig. 1 1 2I >. 
[NGC 2005 is an old globular that has been cited by other 
autho rs as a core-co llapsed object; e.g., Mateo ( 1987p. See 
Mack ev & Gilmorel ( 2003a) for discussion of other possible 
core-collapse candidates in this sample of LMC GCs.] The 
other main outlier in this plot, with /iy,o(W)-/iyo(K) = -1 .09, 
is R 136 = 30 Doradus in the LMC. This very young ob- 
ject's extremely compact configuration is not particularly well 
fit by any of our m odels (for further discussion, see, e.g., 
Mac kev & Gilmorel2003ah . 

The upper middle panel shows a tendency for the projected 
core radius to be larger in the Wilson models for these clus- 
ters, by ~ 5% on average but up to ~ 15% when Wilson 
spheres fit very much better than King models. This might 
be viewed as a demonstration that fitting the "wrong" type of 
model to a cluster introduces possible error of this order in 
the fitted R c . Note that NGC 2005 and R 136 are again out- 
liers here, with their much smaller core radii in the Wilson fits 
corresponding to their much brighter /iy.o- 

It can be seen in the upper right-hand panel of Fig. [^] 
that, when Wilson models fit better than King in the 
LMC/SMC/Fornax sample, the model half-light radius R/, is 
essentially the same — again to within about 5% on average — 
in either model. The reason is that observations out to R > 5Rh 
are generally required ("Fig.1141 in order to show a clear prefer- 
ence for one model or the other. In such cases R/, is very well 
constrained by the data themselves and must be reproduced 
by essentially any model fit. When Wilson- and King-model 
fits have comparable x 1 , on the other hand, it occasionally 
happens that their R^ values differ by factors of 2.5-3. When 
this occurs, it is most often because R\ ist /Ri, is of order 1 or 
smaller, meaning that few if any data are available at large 
cluster radii to constrain the very different extrapolations of 
the two types of model. It then becomes unclear which Ri, is 
correct — although Fig.El s h° ws that our estimated errorbars 
in these cases are also larger than average, properly signaling 
the problem. 

The bottom panels of Fig.^]are all very similar. The pre- 
dicted central velocity dispersions for these clusters are (ex- 
cept for NGC 2005 and R 136) the same to within ~ 1%, on 
average, whether Wilson or King models are fit (and regard- 
less of which is the better fit); total cluster luminosities are 
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FIG. 15. — Comparison of physical cluster properties derived from Wilson-model fits to LMC, SMC, and Fornax clusters, vs. those derived from King-model 
fits. The point enclosed in a square in every panel is the LMC globular cluster NGC 2005. See text for details. 



about as well determined as the cluster half-light radii (with 
the same potential for some large Wilson-King discrepancies 
and large errorbars when the two models have comparable x 2 
values); and global cluster binding energies differ by < 5% on 
average between the two models. 



Figure^]is analogous to Fig.^]but compares the fits of 
power laws to those of Wilson models for the MG03 cluster 
set. NGC 2005 is again enclosed by an open square in every 
panel, and the rightmost datapoint is R 136 = 30 Dor. Sim- 
ilar comments apply to these comparisons as to the previous 



28 



McLaughlin & van der Marel 




I 1 1 1 I 
•— 



I 1 1 1 I 




_L 



_L 



o 



CD 
6 



6 



6 



01 -sf CD I 

odd 

l l l 



[(M) q H/(ld) M H] 3o t 



[U) q 3/(ld) q 3] Sot 





co 
d 



CD 

d 



o 



d 



[(M) ,0, l/(ld) w l] Soi 



U) 0A ^-(id) 0A 7/ 



[(M) 0< V(ld)° 3o T 



X 

+ 



_ OJCL. 



NO. 

X 



- OJ 



CM ^ 

X 

+ 



X 

I 

J 
cm a. 
X 



d 




FIG. 16. — Comparison of physical cluster properties derived from power-law model fits to LMC, SMC, and Fornax clusters, vs. those derived from Wilson- 
model fits. The point enclosed in a square in every panel is the LMC globular cluster NGC 2005. See text for details. 



ones, with the important exception that the half-light radii im- 
plied by power-law fits can sometimes be orders of magnitude 
larger than the R/, obtained from Wilson- or King-model fits. 
This reflects a fundamental difficulty with power-law models: 
extrapolation of a fit which happens to fall too near 7 = 3 over 



some available (too small) range of data is barely convergent 
(see eq. 0) and clearly unphysical. Nevertheless, our esti- 
mated errorbars on R/, (and on L tot ) even for unrealistic power- 
law fits such as these do reasonably reflect the situation. And 
in many of the cases seen here, the power-law fit is, after all, 
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FIG. 17. — Comparison of physical cluster properties derived from Wilson-model fits to Galactic globular clusters, vs. those derived from King-model fits. 
The point enclosed in a square in every panel is Palomar 1. See text for details. 



measurably worse (A > 0) than the best Wilson-model fit. 

Figure flTl next compares the King and Wilson fit param- 
eters for Galactic globular clusters. The point enclosed in 
an open square in every panel is Palomar 1, which has 
Mv,o(W)-Mv,o(K) = 5.81 and log[R c (W)/R c (K)] = 1.74 as a 



result of the much higher concentration of the King model 
fit in this case (see Fig. ^] and recall the discussion around 
Fig. ^] above). The other main outliers in the upper left- 
hand panel of Fig. [17] are NGC 5272 = M 3, at A = -0.87 
and /iv : o(W)-/iv,o(K) = -0.51, and Palomar 2, at A = -0.23 
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and //y,o(W) - /^y.o(K) = -0.72 with a large errorbar. Aside 
from these objects, Wilson fits to Galactic GCs tend to imply 
central surface brightnesses slightly fainter than King-model 
fits, but by only 0.02 mag on average and less than ~ 0.1 mag 
in most cases. Wilson models are also associated with slightly 
larger projected half-intensity radii R c and smaller half-light 
radii R/„ and thus relatively lower concentrations as measured 
by the ratio Rh/R c - This is simply a result of the more ex- 
tended cluster structure essentially assumed ab initio in the 
model. It is apparent again that the disagreement between 
King- and Wilson-model estimates of Rh, and subsequently 
L tot , is potentially largest when the two models fit a given 
cluster about equally well. As above, equal-quality fits with 
disparate Rh typically occur when the available cluster data 
do not extend far enough in clustercentric radius to defini- 
tively constrain the large-/? extrapolation of either model; but 
our estimated uncertainties on the derived quantities gener- 
ally reflect this. The predicted central velocity dispersion and 
global cluster binding energy are as well-behaved as in our 
LMC/SMC/Fornax cluster fits. 

Finally, in Fig. ^] we compare the extrapolated tidal radii 
from King- and Wilson-model fits to 37 of the LMC, SMC, 
and Fornax clusters (top panel) and 43 Milky Way globu- 
lars (bottom panel). We have only included clusters in these 
graphs if they have R\ as t/Rh > 4 — so that the extrapolations 
of the fits to r, are constrained as well as possible — and 
A = (Xwiknn - xk g )/(Xwii.nn +xk g ) < 0.1— in which case 
Wilsorj (119751) s pheres fit the surface-brightness data at least 
as well as Kins ( 1966) models. 

What we have actually plotted in Fig. ED is the mass- 
normalized tidal radius evaluated for each of these clusters 

—1 /3 

within each of the two model fits: r,M t J , the inverse cube 
root of the average cluster density. This quantity is related di- 
rectly to the tidal field in which a cluster is embedded; most 
simply, in the case of a cluster at radius r gc in a spherical 
galaxy, M tot /r f 3 oc Mg a i(r gc )/r gc . Ideally, we would like to 
—1/3 

compare estimates of r,M tot from fits of structural models, 
to the value expected for any cluster in a given galactic tidal 
field (in order, for example, to assess whether an extended 
cluster halo is "unbound;" cf. SI4.3.2> . However, to do this 
requires detailed knowledge of the parent galaxy mass pro- 
file M(r gc ) (including dark matter); of the instantaneous three- 
dimensional position of the cluster, r ec ; and of the shape and 
energy of the cluster orbit, which together set the coefficient 
connecting the mea n cluster density to the average g alaxy 
densi ty at r gc (e.g., iKingl (19621: llnnanen. Harris. & Webbinkl 
119831) . Estimating all of these quantities individually for each 
object in our sample is clearly out of the question, and even 
a statistical treatment of the cluster ensemble in each galaxy 
(such as in llnnanen. Harris. & Webbinkll 19831) is beyond the 
scope of our analysis. 

We simply point out that — as expected — the mass- 
normalized tidal radii implied by the Wilson-model fits to our 
clusters are systematically larger than those implied by King- 
model fits. The dotted line in each panel of Fig.ll8lindicates 

—1 /3 

equality between the Wilson and King values for r,M tot , 
while the bolder, dash-dot lines show the median ratios of 
the two estimates: ~ 2.9 in the LMC/SMC/Fornax sample, 
and ~ 2.5 for the Milky Way globular clusters. We emphasize 
again that each of the points plotted represents a cluster which 
is fit at least as well or better by a Wilson sphere vs. a King 
model. But the question remains open as to whether the lim- 
iting radii (or mean densities) of the former models are quan- 
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FIG. 18. — Mass-normalized tidal radii inferred from Wilson-model 
fits, vs. those from King-model fits, for 37 clusters in the LMC, SMC, 
and Fornax (top panel) and 43 globular clusters in the Milky Way (bottom 
panel). Only clusters with R Vi!t /R h > 4 and A = (x wilson - X King )/(X W ii S on + 
Xicing) ^0.1 Me plotted. Dotted line in each panel indicates equality, 



r r M 10 ' t '' 3 (Wilson) = r,A/ to l / 3 (King). Bolder, dash-dot lines are at the median 
ratios /•,M to 1 / 3 (Wilson)//-,M to '/ 3 (King) = 2.9 for the LMC/SMC/Fornax sam- 
ple, and r t M^ t |/3 (Wilson)/r f M to 1 / 3 (King) = 2.5 for the Milky Way globulars. 



titatively consistent with a naive association of fitted r, values 
with the true tidal radii of the real clusters. The large errorbars 

—1 /3 

on both King and Wilson values of r,M tot ' for all the clusters 
in Fig.^Donly stress further the difficulty of precision work 
along these lines. Even for these best-observed clusters, r t 
is almost always inferred from extrapolation rather than mea- 
sured directly; it is the most uncertain cluster parameter that 
we estimate. 

5. OBSERVED VELOCITY DISPERSIONS AND DYNAMICAL 
MASS-TO-LIGHT RATIOS 

As we have described, the derivation of total cluster masses 
(and all the dependent quantities in Table [Q] above) from 
our surface-brightness fits has been facilitated by adopting a 
mass-to-light ratio for each cluster based on the population- 
synthesis modeling of This is a necessary step towards 
examining physical trends and dependences among the prop- 
erties of star clusters spanning a wide range of ages. It re- 
mains to be shown, however, that the population-synthesis 
models we use predict mass-to-light ratios that are consistent 
with what can be inferred directly for the minority of clusters 
which have measured velocity dispersions. Such a demonstra- 
tion is the purpose of this Section. 

To make this check, we have compiled velocity-dispersion 
data for as many of our modeled LMC, SMC, Fornax, and 
Milky Way clusters as we could easily find in the published 
literature. We have not made an attempt at a comprehensive 
collection, but simply one including enough clusters to ad- 
dress meaningfully the question at hand (in fact, for Milky 
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Way globulars we hav e relied exclusively on the work of 
iPrvor fe Mevlanl lT993l which is itself a compilation of ear- 
lier studies). A drawback is that very few young Magel- 
lanic Cloud clusters have well-determined velocity disper- 
sions: among those in our sample, we have found data only 
for NGC 1850, 1866, 2157, 2164, and 2214 in the LMC, 
and NGC 330 in the SMC. Although bright, these objects are 
usually of relatively low mass compared to the average old 
GC, and their velocity dispersions are intrinsically low and 
difficult to measure. Thus, our comparison of population- 
synthesis and dynamical mass-to-light ratios really speaks 
most clearly to the old-age limit of the models; but the good 
results in that limit are encouraging. 

Our analysis is detailed in Tablell3l After the cluster name, 
Column (2) of this table repeats the V-band mass-to-light ra- 
tio predicted by population-synthesis models [from Column 
(6) of Table |3 or Column (2) of Table [TH. Column (3) is 
the observed velocity dispersion in the cluster, as reported in 
the literature. Some of these dispersions are based on radial- 
velocity measurements of individual stars spread throughout 
the cluster; others, on integrated-light spectroscopy within a 
finite slit width. Either way, every cr p , b s is in effect a weighted 
average of the cluster's projected velocity-dispersion profile 
over some area on the sky. From the details of each ob- 
servation in the original papers, we have estimated the ef- 
fective radius of a circular aperture with roughly the appro- 
priate area. [Note that this is always R dp = for the Milky 
Way globular clusters, since the observed dispersions in this 
case have already been extrapolated to their central values by 
JPrvor & Mevlanl ll9931.1 This is reported in Column (4) of 
Tablell3l Column (5) gives the reference to the source of the 
data. 

Given a King, Wilson, or power-law model with fit- 
ted Wo (or 7) and ro for any cluster, solving Pois- 
son's and Jeans' equations and projecting along the line 
of sight yields a dimensionless velocity-dispersion pro- 
file, <7 P = <r p (R) I <r p (R = 0) as a function of projected clus- 
tercentric radius R. The weighted average 6> 2 (f? ap ) = 

RI(R) a\ dR ap RI(R) dR ' then gives the predicted 

mean-square velocity dispersion within any circular aperture 
of radius 7? ap . We have calculated S for each of the model 
fits to each of the clusters in Table II 3 1 given the aperture radii 
estimated in the table, and obtained the line-of-sight velocity 
dispersions at the cluster centers as (J p {R = 0) = a p , b s /S(R ap ). 
In general, the value of <r p (R = 0) depends on the model used 
to compute S, but the differences in our case are usually small 
and in Column (6) of Tablell3lwe report only the mean of our 
three determinations. 

The observed (J p (R = 0) values are to be compared with the 
predicted er p ,o, based on our population-synthesis M/Ly ra- 
tios for each cluster, in Table 1 121 As described above for the 
calculation of these predictions, our fitted models with known 
Wq or 7 also provide the dimensionless ratio cr p (R = 0)/ (To, for 
(To the theoretical scale velocity appearing in equation (|4}. We 
therefore compute <to and use our fitted ro (TablefTol in equa- 
tion (|4j to compute the central mass density p(r = 0) of every 
cluster in Table ^] A "dynamical" estimate of the V-band 
mass-to-light ratio follows immediately as T v yn = p(r = 0)/ jo, 
with the luminosity density jo taken from Table Our 
estimates of T v ya for every model fit, and the comparisons 

A(log T v ) = log(Ty yn JT p v ° p ), constitute the rest of Table[l3] 
Inspection of Tablell3lshows. first, that the dynamical esti- 
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FIG. 19. — Ratio of dynamical V-band mass-to-light ratio to population- 
synthesis model prediction, as a function of cluster metal licity, for all clusters 
with measured central velocity dispersions in Table ll3l The dynamical Tv 
used are those calculated from cr„ (, s using Wilson-model structural fits to 
each cluster. Population-synthesis mass-to-light ratios are those predicted by 
the model of lBruzual & Charloll 12003ft using the disk-star IMF of lChabriej 
The bold, dash-dot line indicates the median Y^ yn /Tj, 015 = 0.82, 
which has a standard error of ~ ±0.07. 



mates of Tv are generally very model-independent. The sin- 
gle greatest exception is NGC 2005 in the LMC, the cluster 
for which our Wilson-model fit has a much higher central sur- 
face brightness and smaller scale radius, and thus a smaller in- 
ferred mass-to-light ratio, than the King- or power-law model 
fits (see Fig.ll5land Fig. ll2> . Second, the comparison between 
Y y yn and T v op is favorable. Figure[l9]shows this graphically, 
with the ratio of the two mass-to-light values plotted as a func- 
tion of cluster metallicity (from Table [8] above). For definite- 
ness, the results for T v yn from our Wilson modeling have been 
used in this plot, but it makes no significant difference if the 
King- or power-law model numbers are used instead. In all 
cases, the median Ty^/T^ ~ 0.82 ± 0.07. 

The few young massive clusters in the LMC and SMC for 
which we have obtained dynamical mass-to-light estimates 
are shown as open symbols in Fig. \19\ Although these tend 
to fall nominally above the line T v yn = T v ° p , their measured 
c p ,obs m Tablell3lare relatively uncertain, as our associated er- 
rorbars attest. In fact, in three of the six cases only upper lim- 
its to cr p ,obs are given by the original authors. All other points 
in Fig.ll9lrefer to old (r > 10 10 yr) globular-type clusters and, 
as mentioned above, provide a direct check only on that ex- 
treme of the population-synthesis models. Overall, however, 
we feel confident that our use of T v ° p in general to infer mass- 
based cluster properties from simple surface-brightness mod- 
eling is well justified. 

6. k-SPACE PARAMETERS AND GALACTOCENTRIC DISTANCES 



We anticipate a main use of our results in Tables FTOlll II and 
above to be in the definition and interpretation of correla- 
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tions between the primary physical properties of star clusters. 
Ultimately, such correlations can constrain theories of cluster 
formation and evolution. They have been identified and dis- 
cussed in many forms in the literature for old, globular clus- 
ters in the Milky Way and a few other galaxies. So far as we 
are aware, our work here is the first to allow for systematic 
investigatio n of the effec ts of fitting GCs with models other 
than that of lKind |[l966). It is also the first to put a signifi- 
cant number of young massive clusters on a completely equal 
footing with the old globulars. 

Correlations among GCs are typically couched in terms of 
a structural "fundamental plane" a nalogo us to that orig inally 
define d by iDio rgovski & Davisl lll987l) and iDressler et alJ 
( 119871) for dynamically hot galaxies. There are at least three 
equivalent formulations of the globular cluster fundamental 
plane in the literatur e: 

First, iDiorgovskH (119951) presents strong bivariate correla- 
tions involving <t Pj o, ro, A*v,o, Rh, and showing Galactic 
GCs to be an essentially two-p arameter family. 

Second, McLaughlin (2000) works with Ty, L tot , King- 
model c, and global binding energy Ef, to arrive, in different 
form and with different physical emphasis, at the same ba- 
sic conclusion. McLaughlin shows explicitly the equivalence 
between his and Djorgovski's formulations of the GC funda- 
mental plane. In either, the Galactocentric positions R sc of 
the globulars are an important external influence; it is well 
known, for example, that GC half-mass radii and binding en- 
ergies correlate significantly with t heir location in the Galaxy 
(Rh oc R°„r and Eh oc /C° 4 : |van den Bergh. Morbev. & Pazder 
fl99ll lMcLaughlin 2000). 

To bring young massive clusters, such as those we 
have modeled in the LMC and SMC, into analyses along 
these lines, it is preferable to work in terms of M tot , Eo, 
and £/, — rather than their luminosity or surface-brightness 
equivalents — so as to avoid purely age-related effects. All but 
one of the required fundamental-plane variables for our full 
cluster sample are then given in Table s ITOl throu gh 1 1 2 1 above . 
The last — cluster positions within their parent galaxies — is 
listed in Table 1 141 discussed below. 

A thi rd equivalent f ormulation o f the fu nda mental plane is 
that of iBender. Burstein. & Faberl ill 9921) and Burstei rTet al.l 
( 119971) . who manipulate the basic observables of velocity dis- 
persion, surface density, and half-mass radius to define an or- 
thonormal set of derived parameters, 

K\ = (\oga^ fi + logR h )/V2 

k 2 = (log <72 + 2iogE A -log R h )/ V6 (5) 

k 3 = (loga* -log'E h -logR h )/V3 

and find tight distributions in kj, vs. K\ for early- 
type galaxies and (separately ) fo r globular c lusters. 
IBender. Burstein. & Faberl ( 119921) and iBurstein et all (11991) 
actually define this "k space" using the luminosity intensity 
Ih averaged over the half-light radius R^; but in order to 
remove the influence of age from comparisons of clus- 
ter structures, we instead use the average mass density 
Sft = T v h = M tot /2nRf r Then K\ «-> log (cr? „/?/,) is related 

to the total mass of a system, and K3 <-* log(<jp Ri, / M tot ) 
contains the exact details of this relationship — that is, infor- 
mation on the internal density profile. In fact, the mass-based 
«3 of equation (0 can be viewed as a replacement for King- 
or Wilson-model concentrations c or power-law indices 7, 
or any other model-specific shape parameter. As such, any 
trends involving K3 are directly of relevance to questions 



concerning cluster (non)homology. The definition of k 2 
is chosen simply to make the three k parameters mutually 
orth ogonal; it results in the corres pondence k 2 <-> log(X^) 
(see Bender, Burstein, & Faber 1992, for further discussion). 

Tablell4lgives the values of K\, k 2 , and K3 for all the y oung 
and globular clusters that we fit with structural models in 34.21 
In calculating these parameters, we have used the <r p ,0 and £/, 
values as predicted in Tablell2lbv our adoption of population- 
synthesis mass-to-light ratios. Equations (0 are evaluated for 
(Tp O in units of km s" 1 , £/, in M pc" 2 , and Rf, in kpc (fol- 
lowing Bender et al. and Burstein et al., who originally had 
galaxies in mind). 

Table [2] also contains the observed distance, in kpc, of 
each cluster from the center of its parent galaxy. These are 
projected galactocentric radii for the LMC, SMC, and Fornax 
clusters, and three-dimensional radii for the Milky Way glob- 
ular clusters. For the Galactic GCs , we have sim ply copied 
R gc directly from the catalogue of iHarrisl (119961) . For the 
LMC, SMC, and Fornax clusters, we have computed R gc our- 
selves, using th e right ascensions an d declinations of the clus- 
ters as given by Mackev & Gilmore ( 2003a b c) and taking the 
galaxy centers to be 

a 2 ooo = 5 h 25 m 6 s S 2m = -69° Al 1 (LMC) 
a 20 oo = h 52 m 45 s 5 2m = -72°49'43" (SMC) 
a 2 ooo = 2 h 39 m 59?3 <5 200 o = -34°26'57" (Fornax) . 

(6) 

The HI and optical centers of the LMC are offset from each 
other, and here we have mor e or less arbitrarily a dopted the 
optical center of the barfromjvan der Marel ( 2001 ). The SMC 
center is taken from Westerlund ( 1997), and the Fornax center 
from SIMBAD. In converting angular distances to kpc, we 
assume a constant distance of 50.1 kpc to all clusters in the 
LMC; 60.0 kpc to the SMC; and 137 kpc to Fornax. 

7. SUMMARY 

We have fit three distinct dynamical models to V-band 
surface-brightness profiles for each of 68 massive star clus- 
ters (50 of which have young ages, between several Myr and 
a few Gyr) in the LMC, SMC, and Fornax dwarf spheroidal, 
and to 85 old globular clusters in the Milky Way. We have 
also applied publicly available population-synthesis models 
to infer the expected intrinsic (B-V)o colors (and thus red- 
denings) and V-band mass-to-light ratios for all of these clus- 
ters plus another 63 Galactic globulars. Combining these with 
the surface-brightness model fits, we have calculated a wide 
range of structural and dynamical parameters characterizing 
the clusters. Our main results are contained in Tables 151 1 101 
^2E1 and^^above. We have also taken velocity-dispersion 
measurements from the literature for a subset of the full clus- 
ter sample and calculated dynamical mass-to-light ratios for 
comparison with the predicted pop ulation-synthesis values. 
The results of this are in Table 1131 which shows quite good 
agreement in general. 

The V-band surface-brightness data we have employed 
for LMC, SMC, and Fornax star clusters derive primarily 
from H ST-base d starcounts m ade in the inner regions by 
Mackev & Gil morel ( l2003albld collectively MG03). We have 
supplemented these wherever possible with ground-based 
starcounts and BV aperture photometry from the literature, 
both to extend the MG03 cluster profiles to larger projected 
radii and to re-calibrate the V-band magnitude scale of MG03. 
This re-calibration turns out to be rather significant, amount- 
ing to several tenths of a magnitude (in the sense that the sur- 
face brightnesses published in MG03 are generally too faint) 
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in many cases. Full details of our analyses on the 68 MG03 
clusters are in 32. 11 the re-calibrated surface-brightness pro- 
files, with all ground-based data properly incorporated, are 
given in Table [5] All surface-brightness profiles for Milky 
Way globular clust ers were taken from the database con- 
structed by Trager. King. & Dio rgov ski ( 1995), with only mi- 
nor modifications described in 32.21 

Our population-synthesis modeling, wh ich includes a com- 
parison of results from two separate cod es ( Bruzual & Chariot 
120031 iFioc & Rocca-Volmerangell9 97) using a variety of as- 
sumed stellar IMFs, is described in |5] Table [8] lists intrin- 
sic cluster colors and theoretical mass-to-light ratios obtained 
with six different code+IMF combinations, although for our 
subsequent modeling (requiring estimates of cluster V-band 
extinctions and conversions betwee n luminosity and mass) we 
adopted numbers from the codeofj Bruzual & Chariot (2003) 
using the disk-star IMF of Chabrierl ( 120031) . As described in 
33.21 we have thus obtained systematically lower predicted 
mass-to-light ratios, at any age, than MG03 inferred for their 
LMC/SMC/Fornax clusters. Our values nevertheless compare 
well with the dynamical mass-to-light ratios calculated di- 
rectly from velocity-dispersion measurements for 57 clusters 
(mostly old globulars) in |5] Additionally, our model-based 
extinctions agree very well with direct measurements for the 
old globular clusters in the Galaxy and the Fornax dwarf. 

The three models that we fit to each cluster are described 
in some detail in which also contains the fits themselves 
and the bulk of our derived structural and dynamical param- 
eters. The models are: (1) the s i ngle-m ass, isotropic, mod- 
ified isothermal sphere of iKing l {1966); (2) an asymptotic 
power law with a constant-density core; and (3) an alternate 
modification of the isothermal sphere (still single-mass and 
is otropic) based on the stellar distribution function developed 
bv lWilsonl d 1 9751) . For otherwise similar clusters (e.g., given a 
fixed total lumin osity, c entral surface brightness, and effective 
radi us). |Wilsonl i 19751) spheres are spatially more extended 
than King ( 1966) models, although both have finite tidal radii; 
the untruncated power-law models we fit are formally infinite 
in extent and, in some cases, barely convergent in their inte- 
grated properties. Even so, the structural differences between 
these models are most important at relatively large cluster- 
centric distances, in the outer halos beyo nd a few effective 
radii. As a result, we have verified ( 34.3> that for most clus- 
ters, most of the physical properties we have calculated are 
reasonably well constrained no matter which particular model 
is taken to fit the surface-brightness data. 

We have looked with particular interest at the question of 
which of the three fits — the spatially limited King model; 
a more extended but still finite Wilson sphere; or an infi- 
nite power-law model — pro vides the best description of these 
clusters in a x 2 sense (see 34311 . All of the 68 LMC, SMC, 
and Fornax clusters from MG03 are fit at least as well, and 
in several cases significantly better, by the larger env elopes 
oflWilsonl ill 9751) models rather than the more compact lKing l 
(1966) models. This is true for the 18 old globular clusters 



as well as of the 50 younger objects in this sample. It also 
holds for the majority (w 70/85) of the old Galactic globulars 
to which we have fit all three mo dels; a particularly clear ex- 
ample is provided by ui Centauri ( 34.2. 11 1. In all cases, asymp- 
totic power laws are not vast improvements over the Wilson- 
model fits; in fact, they are somewhat worse for many clusters. 

In the LMC, SMC, and Fornax sample especially (where 
the surface-brightness data are most homogeneous), there is 
no correlation between cluster age and the relative quality of 
fit for King vs. Wilson or power-law models. Instead, we 
have shown that the primary factor in determining whether 
an extended-halo model describes a cluster better than a King 
(1966) model is simply the spatial extent of the available 
surface-brightness data being fit. Specifically, it is only when 
a cluster's observed density profile reaches to more than ~ 4— 
5 effective radii that it becomes possible to decide conclu- 
sively whether or not it has a non-King envelope structure; 
and when such data do exist, a more distended but finite Wil- 
son model is most often the better option of the three we have 
examined — whatever the cluster age. 

Thus, we conclude that the extended halos which are known 
to surround many young massive star clusters in the Mag- 
ellanic Clouds and other galaxies do not require description 
by untruncated power laws which would necessarily have all 
the clusters ove rfill the Roche lobes defined by their par- 
ent galaxies (see lElson. Fall. & Freemanl! 19871 lLarsenll2004t 
Schweizer 2004) — although whether or not this does happen 
in individual cases is a complicated question that we have not 
undertaken to address here. More generally, despite the ad hoc 
nature of the Wilson and power-law models that we have fit, 
it is a clear fact that self-gravitating clusters commonly have 
envelo pe structures which do not match the extrapolations of 
simple King ( 1966) models fitting the cluster cores. This phe- 
nomenon is not confined exclusively to young clusters and is 
not obviously only transient; it may point instead to generic, 
internal cluster physics not captured by King's stellar distri- 
bution function. 

This interesting issue aside, the work we have presented 
provides the basic information required to define and compare 
physical parameter correlations and the fundamental plane(s) 
of young and old massive star clusters in the Milky Way, 
LMC, SMC, and Fornax. It also offers a starting point for 
careful examination of potential model-dependent artifacts in 
these correlations, and ultimately it should allow for direct 
contact to be made with the well-established fundamental 
plane of elliptical galaxies and bulges. We plan to address 
these issues in future work. 
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ONLINE MATERIAL AND FIGURE 12 

The next seven pages show extracts from Tableland Tables lT0ljl4l which are discussed in 33.21 34.2.21 ^5] and 50 Each of 
these tables can be downloaded in its entirety, in machine-readable format, from the online edition of the Astrophysical Journal 
Supplement Series, where this paper is published as 

McLaughlin, D. E., & van der Marel, R. P. 2005, ApJS, 161, 304. 

Following the tables is Figurefl2l( 34.2.2> . which shows the full set of fits to each of the LMC, SMC, Fornax, and Milky Way 
star clusters that we have modeled in detail. The full-resolution versions of these plots are available on the ApJ website, and also 
from http://www.astro.le.ac.uk/ ^dml31/clusters.html . 



TABLE 8 

Population-Synthesis Colors and Mass-to-Light Ratios 



Published Data a 



SSP Model Code/IMF b 



^.luster 


age 


[rc/rij 


\tS—V ) 


Property 


Ijv /v Illy 




ppr./rhn 


x^IHj/Ci<Aip 




rCj\Ji JS. 1 kj 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


LMC-HODGE11 


10.11 


-2.06 


0.64 


(fi-V)o 


0.594 ±0.015 


0.616±0.014 


0.609 ±0.017 


0.626 ±0.016 


0.606 ±0.017 


0.610±0.017 










1 V 


1.893±0.160 


3. 133 ±0.226 


2.013 ±0.166 


3.279 ±0.237 


1.778 ±0.126 


1.923 ±0.146 


SMC-KRON3 


9.78 


-1.16 


0.68 


(B-V)o 


0.664 ±0.017 


0.675 ±0.017 


0.700 ±0.021 


0.706 ±0.021 


0.699 ±0.020 


0.698 ±0.021 










^pop 


1.129±0.113 


2.018±0.194 


1.281 ±0.140 


2.269 ± 0.240 


1.236±0.130 


1.312±0.127 


FORNAX 1 


10.11 


-2.20 


0.65 


(B-V)o 


0.579 ±0.015 


0.601 ±0.015 


0.589 ±0.017 


0.608 ±0.016 


0.586±0.017 


0.590 ±0.017 










-y^pop 


1.915 ±0.163 


3.155±0.228 


2.004 ±0.167 


3.249 ±0.239 


1.763 ±0.127 


1.911 ±0.147 


MW- 1636-283 


10.11 


-1.50 




(B-V)o 


0.659 ±0.016 


0.678 ±0.016 


0.698 ±0.025 


0.710 ±0.024 


0.694 ±0.025 


0.698 ±0.025 










^pop 


1.888±0.155 


3. 169 ±0.224 


2.134±0.187 


3.542 ±0.286 


1.914±0.149 


2.058 ±0.170 



NOTE. — A machine -readable version of TableHjlis published in its entirety in the electronic edition of the Astrophysical Journal Supplement Series; only a portion is shown here, for 
guidance regarding its form and content. 

a Column (2) is th e logarithm of cluster age in years. Ages, metallicities, and their uncertainties for LMC, SMC, and Fornax clusters are taken from the compilation and analysis by 
Mackev & Gilmore i 2003a be). All ages given as greater than 13 Gyr by Mackey & Gilmore have been re-set here to 13 ± 2 Gyr (log r = 10. 1 1 ± 0.07). Aperture (B — V) colors for 
LMC, SMC, and Fornax cluste rs are from the literature summarized in Table|3l Milky Way globular cluster ages a re set to a unif orm 13 ± 2 Gyr (log r = 10.1 1 ± 0.07); their metallicities 
are taken from the catalogue of Harris < 1996) and assigned uncertainties of ±0.2 dex. (B — V) colors are also from Harris < 1996), where available. 

Key to the combinations of population-synthe sis codes and assumed stellar IMFs: BC/ChD — code of Bruzual & Chariot < 2003) using the disk IMF of Chabriei i 2003); BC/Salp— 
code of iBruzual & Charloi 120031 u sini^ the ISalpeteil I1955l> IMF; PEG/ChD — PEGASE v. 2.0 code of lRoT^^occ^^lmer angel 119971) using the disk IMF of lChabrieil 1200311 : 
PEG/Salp— PEGA SE code using the |Salpete!l Il955l> IMF: PEG/ChGC— PEGASE code using the Globular Cluster IMF of |Chabriei 2003); PEG/KTG— PEGASE code using the 
IMF o flKroupa. T~ & Gilmore ■ 1993T 
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TABLE 10 
Basic Model Fit Parameters 



Cluster 




A v 


D [kpc] 


Wpts 


weight 


model 


9 

X 


W / 7 


c 


/iyo 


ro [sec] 


ro [pc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


LMC-HODGE 1 1 


"•"-0.05 


U - -0.08 


50.1 


50 


1 


K 


17.62 


fS 


1 ' 2:, -o;i4 




14.42 + ^''^ 
l5.2l + jS 


i SO 














W 


17.06 


f\ 1 +0.8 




19.32+0; ' 


-1 sr\-\-i) V> 

3 - 69 ^8:l5 


SMC-KRON3 


u ' 1J -0.05 


05 +0 08 
u - UJ -0.08 


60.0 


72 


1 


PL 
K 


16.73 
56.99 


3 65+ 35 
5 5^ 


00 

j . 4+0.05 


19.32+010 
20 04+^°' 


16.52^:40 
23.24+j : p 
28.15_j'c,j 


4 01+0-58 

6 76+^ 

"■'"-0.38 














W 


67.83 


4 4+u.j 
^-0.2 


^'-0.05 


zu - uy -o.io 


8 1 9+0 -38 














PL 


87.98 


5 10+ 025 

-oas 

2 6+°-5 
z -°-0.6 


OO 


20 13+ 009 


37.46+f™ 
21.50?| 


! 0.90+0 .so 
14.28+f™ 


FORNAX 1 





a 99+O.O8 
u ' z -0.08 


137.0 


51 


1 


K 


31.25 


61+° 08 
uol -0.10 


2S ofi +0. ID 














W 


32.26 


i+°- 7 


a, +0.49 

u - ul -o.oo 


23 07 +0.08 


116.30^ 


77.24 + o i ,°o 














PL 


33.12 


q 7 c+2.45 


OO 


23 1 5+ ' 09 


39.6l4f? 


26 31+5 05 
ZD ' -3.46 


MW-AM1 








121.9 


43 


1 


K 


44.00 


6 6+ 04 


1 41+OI2 


3 S O+0.04 
/J - SS -0.03 


lz - 1I -0.82 


7 I6+ - 69 
'• lo -0.49 














W 
PL 


45.52 
46.16 


6 6+° 6 

-0,2 

3 451+0.36 
-0.20 


9 1 1+0.43 
z " 11 -0.36 
OO 


93 QO+0 °4 


12.92+ 147 
13 43+-^ 


7.64+0.87 

7 94+^:^ 
'• 7 ^-0.98 



NOTE. — A machine -readable version of Table l 1 Ol is published in its entirety in the electronic edition of the Astrophysical Journal Supplement Series; only a portion is shown 
here, for guidance regarding its form and content. Key to columns: 
Column (1) — Cluster name. 

Column (2) — Zeropoint correction applied to published V-band surface brightnesses; see jKlfor details. 

Column (3) — V-band extinction, from population - synthesis modeling for LMC and SMC clusters (Table with Ay = enforced for clusters with observed (B — V) bluer than 
model (5 — V)o); from Mackev & Gilmore 2003c) for Fornax globular clusters; from Harris 1996) for Galactic GCs. 
Column (4) — Heliocentric distance 

Column (5) — Number of surface-brightness datapoints constraining model fits 

Column (6) — Weighting scheme for model fitting: "I" denotes weighted least-squares fit to surface -brightness data in linear units (Lq pc~ 2 ); "s" denotes weighted least-squares fit 
to surface brightness data in mag arcsec -2 ; "u" denotes unweighted least-squares fit to surface brightness data in mag arc sec -2 . 

Column (7) — Identification of fitted model: K- Kins 1966 ) model, defined by equation W=isotropic Wilson 1975) model defined by equation j^J; PL=power-law with core, 
defined by equation j2|. 

Column (8) — Minimum x 2 f° r ° es t model fit. For weighting schemes '1" and "s," this is the unreduced x 2 °f me best-fitting model; the value per degree of freedom is given by 
X 2 /(Afpi s -2) with 7Vpt s in Column (5). For unweighted fits "u," x 2 * s actually the rms deviation of observed surface brightness (in V mag arcsec - ") from the model profile. 
Column (9) — For King- and Wilson-model fits, entry is the dimensionless central potential Wo; for power-law models, entry is the exponent 7 of the asymptotic three — dimensional 
density profile (see eq. |2*1). 

Column (10) — Concentration parameter c = log(r, / ro) (undefined for power-law models, which have infinite spatial extent). 

Column (11) — Best-fit central V-band surface brightness, after correction for any zeropoint change in Column (2) and extinction in Column (3). 

Column (12) — Best-fit model scale radius ro in arcsec. 

Column (13) — Best-fit model scale radius ro in pc. 



TABLE 11 
Derived Structural Parameters 



Cluster 


mod 


log n 


log R c 


logR;, 


\og(R,,/R c ) 


log ib 


log 70 


log Liot 


Vm 


log 








[pc] 


[pc] 


[pc] 




[i Q pc- 2 ] 


[i pc- 3 ] 


[Z-0] 




[L pc- 2 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



LMC-HODGE1 1 



SMC-KRON3 



FORNAX 1 



MW-AM1 



K 
W 
PL 
K 
W 
PL 
K 
W 
PL 
K 
W 
PL 



1.80+ 



n+0.13 
u -0.11 
? 41+0-43 

OO 
+0.02 
-0.02 



1.97 
2.18+ 



"-0.03 
OO 

1.76+0- 03 

, on+Sffi 

1 - oy -o.oo 
oo 

9 97+0.08 

-QP8 
9 99+0.39 

OO 



510+O-023 
U - Jiu -0.021 
n s I 8+0.020 
U - MS -0.023 

522+0- 025 

-fl.flja 

n 707+0.019 
-8.021 

0- 81721 
0.840+00.6 
0.999+0 026 

1 oog+o^Oti 

1 n37+0019 

1 - UJ '-0.0I8 

0.8294*j 
0.843|»| 
0.840+00 



0.843+0063 

900+°' 135 

1.065+|g 

1 039+0.000 

1 n93+0.001 

1.094 
1 093 

1 - u ' J -0.002 

I nqq+0 00 3 

""-Q.QQ2 
i 245+0.044 
L -*** -0.038 
1 311+0-134 

1.558" 



-0 224 



n aa^+0.083 
U - JJ -0.065 
n 3S1+01 58 
U - iS ^0.083 
.517+0.282 

n 97«+8;o^i 
-0,019 

a 99 9+0.020 
-0.012 

0.183+0-0.3 
0.095+;! 
0.084+fg 
067+0 0' 6 
0.416+o : o™ 
0.468+0;|64 
0.718^ 



9 oc+0.04 
z -° -0.02 

2 84+ 04 
z -°^-0.02 

2 84+0- 04 
-fl-fi? 
9 rc+0.04 

Z - J -0.02 

9 iro+0.04 

J -0.02 
9 C9+0.04 
Z - J -0.02 
1 34+O.04 

-002 
i 34+0.01 
1 - J -0.02 
i n 1+0.04 
1 - J1 -0.02 

1 02+0- 01 
-0Q1 
1 01+ 001 

1.01 



2 03+ 05 

2 0I+0- 05 
z - ul -0.05 

2.01+ 005 

! 46 +8:8§ 

1 40+°' : ^ 

i - HU -0.04 

L 3 6 +o.04 



■II 01 



0.03 

01+ K8i( 
u - ul -0.03 
_ n ()i +0.05 

-0 12+0- 04 

u ' lz -0.05 
_ .14+0.04 

-o 
u -' J -o.06 



4 RQ+0.06 

4 94+0.10 
^■^-0.06 

5 03+0- 13 

-fi-10 
5 10+00+ 

-004 
s ns+0 04 

5 07+ 04 
-OPi 
4.12+003 

4 19+8:83 
-0.03 
4 12+0.03 
^• lz -0.03 
r, 77+O.O2 
J,/ '-O.01 
a 09+0. 06 

3 96+0- 'J 



11 12+ - 12 

1 1 00+0 15 
1 1 .uu^ 24 

10.77+0- 24 

11.00+™ 



M 
Ml 



11.04 
11.07 

15.24; 

1 5 23 

1J ' ZJ -0.08 
1 5 24+0 08 

1 c oc+0.03 

n.«_ 004 

15 72+0.10 

15 39+0 ; « 
,J - J7 -0.38 



2.41+0.06 
2 34+8:89* 

^ ■ 1 U _Q ^ A 

2 17+0 
2.20+0 
2.22+0 



1.13+ 

1.13+0 

1.12+0 

0.48+0 

0.40+0 

0.04+0 



20.40+0- 22 

20.57+o :4 6 

21 03+o ; 5' 
zi.oj „ 4 „ 

21 00+0- 10 
-SOS 
20.91+; 

20.86+; 

23 60+0-08 

23.541 

23.62^| 

25 21 +0 - 18 
95 41+O.51 

-Q.30 

26 3? +2 - 06 



NOTE. — A machine-readable version of Table ll 1 l is published in its entirety in the electronic edition of the Astrophysical Journal Supplement Series; only a portion is shown here, 
for guidance regarding its form and content. Key to columns: 
Column (1) — Cluster name. 
Column (2) — Model fit. 

Column (3) — Logarithm of cluster tidal radius in pc: log r t = c + log jq. By construction, r, = oo for power-law models. 
Column (4) — Logarithm of projected core (half-power) radius at which I(R C ) = /(0)/2, in pc. 

Column (5) — Logarithm of projected half-light (effective radius) in pc. Half of the total cluster luminosity is emitted from within R/,. 
Column (6) — Logarithm of ratio Ri,/R c , a measure of cluster concentration. 
Column (7) — Logarithm of central luminosity surface density in V-band Lq pc~ . 
Column (8) — Logarithm of central luminosity volume density in V-band Lq pc -3 . 
Column (9) — Total V-band luminosity of fitted model. 

Column (10) — Apparent V-band magnitude of model cluster, corrected for extinction: V lol = 4.85-2.5 log (L^/Lq) + 5 /og(D/ 10 kpc). 
Column (11) — Logarithm of V-band luminosity surface density averaged inside the half-light radius: //, = Lt ot /27r^. 

Column (12) — V-band surface brightness corresponding to average projected luminosity density in Column (11): {^v)h = 26.422 — 2.5 log(//,/LQ pc~ 2 ). 
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TABLE 12 
Derived Dynamical Parameters 



Cluster 




mod 


log Aftot 


log£ A 


log S 


logpo 


log £* 




Vesc.O 


log t Th 


log/o 








[M Q ] 


[erg] 


[M pc- 2 ] 


[M oP c- 3 ] 


[M pc" 2 ] 


[kms -1 ] 


[kms -1 ] 


[yr] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



LMC-HODGE 1 1 1.893+° l j«° 1 K 5.17+°<" 49.91+g JJ 2.31+0* 2.69 +0 Q % 3.63+ , 22 14.14+"]' 9.40+ ,-' 2 -0.605+°°« 2 

w c 29 +o:i8 4009+0:11 3 19 +o:o5 2 29+<™ 6 ?6?+° :i ° 3 63+<> :22 14 1 7+0:89 g =0+0:25 _ n 629+0:06* 

vv j.zz_ 007 *y.yt_ ll J.iz_ 005 z.z^_ 00f) z.oz_ |9 j.oj_ 02i i ^- i '_o.82 yJU -0.I2 u ' oz '-0.060 

PL 5.31+°;> 3 49.95!°!' 3.12+°| 2.29+°$ 2.43^ 3.64+°- 22 14.27+^ 9.75+^ -°-642+°;°™ 

SMC-KRON3 1.129+0;JH K 5.15™ 49.65_* : 2.61™ 1.51™ 2.22™ 2.74!; 10.51_g 9.72_g -1.045!«|| 

w c n +0.06 4 o 61+012 0=0+0.06 146+007 2 2fi +0.06 9 78+019 in 42+ 071 9 67+0 03 _ 1n7 +0.056 

-8P.6 -812 -806 -806 -fl.80 -818 -Q-S6 -fl-83 -P-Q41 

PL 5.12!; 49.62!; 2.57!; 1.41™ 2.28!; 2.81!; 10.33!; 9.64™ "1-213™ 

FORNAX 1 1.91S8;li K 4.40« 48.14^;}g 1.63j«;« 0.31™ ^dol l-l^of 4.09™ 9.46™ -1-150™ 

W 4 40+ 005 48 15+ 010 1 62+ 005 29+ 05 I 47+ - 05 1 16+ 07 4 1 +0 23 9 46+ 002 -1 1Q1+0.042 

w ^ u -q.Q5 ^•'-'-O.IO 1 - DZ -0.05 u - zy -0.05 1 -* z -0.05 '-'"-O.Oe ^- lu -0.22 y - w -0.02 1 - lyl -0.024 

PL 4.40+0$ 48 . 14 +o.,o j.59+0.05 . 23 +o.06 L40 +o.05 L15 +o.07 4 . 05 +o.23 9 . 4 7+o.02 -1.266+°|° 

MW-AM1 1.868+°;^ K 4.04™ 47.25™ 1.29™ 0.15™ 0.76*® 0.63™ 2.51+«;jJ 9.55™ -0.473™ 

W 4.10$ : g 47.27-0;8? 1.28_™J 0.13_™ 0.67_™ 0.63!™ 3 2.52+°;j 9.66+^ -0.513!«;"| 

PL 4.23!°;!? 47.29™ 1.28+°;°* 0.14+°$ 0.31^ 0.63+™ 3 2.53!°;! 10.09^" -0.508™ 



NOTE. — A machine-readable version of Table |12| is published in its entirety in the electronic edition of the Astrophysical Journal Supplement Series; only a portion is shown here, 
for guidance regarding its form and content. Key to columns: 
Column (1) — Cluster name. 

Column (2) — Population-synthesis V-band mass-to-light ratio, from Column (6) of Table|8*l 
Column (3) — Model fit. 

Column (4) — Total cluster mass from extrapolation of fitted model density profile. 

Column (5) — Total cluster binding energy. 

Column (6) — Logarithm of central mass surface density. 

Column (7) — Logarithm of central mass volume density. 

Column (8) — Logarithm of mass surface density averaged over the projected half-light/half-mass radius: Sft = M,oi/2tt^, 
Column (9) — Predicted line-of-sight velocity dispersion at cluster center (o"p cc T^. op ). 
Column (10) — Predicted velocity of escape (to infinity) at cluster center (v^ e oc T pop ). 
Column (11) — Median two-body relaxation timescale, evaluated at the half-light radius; see text. 

Column (12) — Central phase-space density fo = /9o/(27rcr?) 3 ' 2 , where a c is the unprojected, one-dimensional velocity dispersion at r = in the model. Units are Mq pc~ 3 km -3 s 3 , in 
which case the core relaxation time is given roughly by log (t ro /yr) C± 8.28 -log fo. 



TABLE 13 

Population vs. Dynamical Mass-to-Light Ratios 



King fits 



Wilson fits 



Power-Law fits 



Cluster 



(1) 



(2) 



^"p.obs 

[kms -1 ] 
(3) 



Rap 
[sec] 
(4) 



Ref. 



(5) 



<t p (R = 0) 
[kms -1 ] 
(6) 



Tf A (log TV) 
[M @ L-^] 

(7) (8) 



ydyn 

(9) 



A(logTy) 
(10) 



-^dyn 

(11) 



A(logT„) 
(12) 



LMC-NGC1466 
LMC-NGC1754 
LMC-NGC1786 



1.88±0.03 
1.88±0.15 
1.87±0.16 



LMC-NGC1835 1.87±0.16 



LMC-NGC1850 
LMC-NGC1866 
LMC-NGC1916 



0.05 ±0.01 
0.10 ±0.02 
1.90±0.16 



LMC-NGC2005 1.88 ±0.16 



LMC-NGC2019 1.87 ±0.16 



LMC-NGC2157 

LMC-NGC2164 

LMC-NGC2210 

LMC-NGC2214 

LMC-NGC2257 

SMC- NGC121 

SMC- NGC330 

FORNAX3 

FORNAX4 

FORNAX5 

MW-NGC104 

MW-NGC1851 

MW-NGC1904 

MW-NGC2419 

MW-NGC2808 

MW-NGC288 

MW-NGC3201 

MW-NGC4147 

MW-NGC4590 

MW-NGC5053 

MW-NGC5139 

MW-NGC5272 

MW-NGC5286 

MW-NGC5466 

MW-NGC5694 



0.05 ±0.01 
0.06 ±0.01 
1.88±0.16 
0.05 ±0.01 

1.87 ±0.15 
1.77 ±0.10 
0.04 ±0.01 
1.88±0.16 
1.72 ±0.08 
1.92±0.16 
2.35 ±0.24 
1.98±0.17 

1.88 ±0.15 
1.90±0.16 
2.02±0.17 
1.97±0.16 
1.88±0.15 
1.87±0.16 
1.89±0.16 
1.93±0.16 
1.87±0.15 
1.88 ±0.15 
1.87±0.16 
1.92±0.16 
1.87±0.16 



3.80 
7.80 
8.30 
9.90 



ft® 



10.40!J;« 
11.50+|| 



3 00+ - 70 
'•""-0.70 

230 -040 

8.20; ' 
16.60 



8.10!!| 
15.30i4;| 



1 I 70+ 3 00 
11 -' U -3.0Q 

7 50+' 30 

—1.30 



o cn+0.60 

4 30+0.00 

7 30 +1 ™ 
3 on+0.00 
-3.90 

5 50+'°° 

-10Q 

5 90+ 2l ° 

-2.-2Q 

8-80-J;?8 

5 10 +1 00 
-1.20 

7 00+ 1 ™ 
'■""-1.70 

1 1 50+ 2 - 30 
11 ••'"-2.30 

10.40;|0| 
+f:83 

3 00 +o: ^ 
-fl-69 
13 40+2 f >« 

2**8 



5.20+ 



,+°:3I 



5.2Q* 
2 60+° 52 
9 cn+0.50 
^'"-0.50 

1.40*| 
16-00!f;|8 
5.60+1J 2 
8.00!J| 
1 70+ ' 34 

'•'"-0.34 

5 50+ 110 
J - -1.10 



2.8 
2.8 
2.8 
1.1 

2.8 
1.1 

40.0 
100.0 
2.8 
1.1 

2.8 
1.1 

1.1 
2.8 

40.0 
15.0 
2.8 
25.0 
40.0 
3.4 
40.0 
1.5 
1.5 
1.5 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 



1 
1 
1 

2 
ave. 

1 

2 
ave. 

3 

4 

1 

2 
ave. 

1 

2 
ave. 

2 

1 

ave. 
2 
5 
1 
5 
2 
1 
6 
7 
7 
7 



3 - 81 -l'il 
7.9 

9.91+ 3 ^' 

^'-o.si 
10.52« : | 

11.52!!;iJ 

3 1Q+0.77 
-'• -0.75 
267+0.47 

i2+?:fl 

6Ji_] 22 

16.64^1 
13.163; 
8.60!|| 

1 5 7 1 +4-44 
1J -' 1 -4.42 

12 67+ 304 
lz - D '-4.06 

I 1 73+ 3 -°l 
7.61+1 33 

9.89+llf 

3 16+™ 

4 55+O.OO 

7 -jo+1.71 
'■ -1.71 

4 1 5+0.01 

5.92+| 7 f 

6 73+0 01 
°-"-6.73 

8.89+J oj 

5 17+ i: " 2 
7.10+ 2 ',* 

II 50 +2:3 ° 
11 "'"-2.30 
in 4(1+2.08 

5 20+ 104 
-104 
3 00+ 0f ' 

13.40+H 
2 91+0:58 

2 c +0.50 
z '-0.50 

1.40*| 

le.oo+J; 20 , 

5.60+jj 2 
8.00!J| 

1 70+ - 34 

'•'"-0.34 
' 10 



1 15+ 2 - 48 
3 .79+5.21 

1.78+81 



2.50 
2.14 



+i:3 



1.42+ 



+?:% 

-1.24 
+1.03 

M 



1+0.50 
-0.41 
■ 71+0.% 

-fl-71 

1 < 6 +o.l4 

12+ 007 
"•' z -0.05 

13+ 007 

"• 1;, -0.05 

0.68 
2.75 
1.72 

2.18_ 073 
7 6I+ 5 - 64 
4.89+ 2 ' 72 
2.97+j 97 

1.25+o : * 
2.11+0:88 

0.22 + ^ : ^ 

6i+8:48 

" D1 -0.61 

2 16+ 1 - 28 
z - io -0.94 

1 00+0 .11 

as 



10.19; 

1.50 
0.54„ 

t 71+I.O8 
'•"-0.96 

, .48+0.80 
4 si* 08 

1 -jo+0.59 
'-"-0.48 

1.61+!!- 7 ' 



1.16; 

0.61 
1.46 
2.15 
2.87 



-1.06 

1.01+ 04 ' 
0.92+ 



m 

^-0.33 
1.1 8+ 053 



2.54+ 



5.5o;; : '° 



°-0.44 
1+1.17 
-0.93 
I 39+0.63 
1 - Jy -0.50 

99+049 
" yv -0.39 

1 61+ 074 
'• ul -0.60 
1 -jc+0.59 
1 - J -'-0.49 



-0.213+0, 506 

304+ - 413 

-0 023+ 8 ' 2 ' 4 
"■""-0.244 

0.05^:1* 

-0 1 1Q+0.170 
"■ lly -0.181 
-0 039+ 0231 

-0 077+ 076 

"■"' ' -0.077 

0.354+O- 2 ' 2 
.ii7+o:»5 

-0.443+^^ 7 ^ 
161* 2 ^ 

-0 044+ ' 243 

0.065+™ 

0.608+gi 7 ! 

0-416+° :23 °, 

202+° :259 
"• zuz -0.314 
_ . 76 +0 .199 
"' '"-0.221 
053+0-' 87 

606 +ii - 2 '° 

0.992+Oj 1 " 

0.059+O; 24 ' 

1.260+« 9 121 

0.736 +g - 201 



mi 



-0.0741 
1.080+]]," 4 
296 +0 - 150 

-0 064+ ' 206 

400+° - 
u.tuu 3U 

-0 947+ 206 
"■ Z -0236 

-0 090+ 8 m 
"■" -0231 

-O2IO+ 8199 
" "-023.1 

-0 497+0.20I 
"■ H; "-0.2.32 

0.038|§ 
-0.269*1 

-0 31S+0-198 

_n 214+0.200 
"■ zl -0.237 

1 32+fl- 202 
o.loz_g 232 

-0.130 _ 228 

-0 275+ - 211 
u - z/J -0.250 

-0 077+ - 203 

-0 142+°- 197 
-0.229 



1.14+2.45 

3.4l41n 



2.53 +: 

2.16; 
1.38 



-0.36 
+0.50 
-0.39 

1 - D -'-0.68 
■ 51+0.14 
1 - 31 -Q.14 

12+ 007 
"•' z -0.05 

12; 007 

0.68*1 

2 71+1.67 

-1-5S 

1 70+ 1 " 2 

1 - lo -0.38 
3.58 +27 « 
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TABLE 13 — Continued 
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King fits 



Wilson fits 



Power-Law fits 



Cluster 



(1) 



^pop 



(2) 



CT p,obs 
[kms- 1 ] 
(3) 



[sec] 
(4) 



Ref. 



(5) 



ff P (R = 0) 
[kms- 1 ] 
(6) 



r dyn 



(7) 



A(logT v ) 
(8) 
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(9) 



A(logTv) 



(10) 



r dyn 



(11) 



A(logTv) 
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MW-NGC5824 
MW-NGC5904 
MW-NGC6093 
MW-NGC6121 
MW-NGC6171 
MW-NGC6205 
MW-NGC6218 
MW-NGC6254 
MW-NGC6341 
MW-NGC6362 
MW-NGC6366 
MW-NGC6388 
MW-NGC6402 
MW-NGC6441 
MW-NGC6535 
MW-NGC6656 
MW-NGC6712 
MW-NGC6715 
MW-NGC6779 
MW-NGC6809 
MW-NGC6864 
MW-NGC6934 
MW-NGC7089 



1.87±0.16 
1.96±0.16 
1.87±0.16 
1.99 ±0.17 
2.09±0.19 
1.88±0.15 
1.89±0.16 
1.88±0.15 
1.93±0.16 
2. 16 ±0.20 
2.28 ±0.22 
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REFERENCES. — (l)-Dubath, Meylan. & Mayo]] fT99^): (2)- Mateo, Welch, & FischeJ f799ll) : ( , 3)- ^sch er. Welch. & Mate j TT991 : C4) 4Fischer etall IT991) : (5) 4Lupton etall fT989l) ; (6) 4Feast & BlacH fT98Cl) ; (7)- 

" JIT 



iDubath. Mevlan. & Mavo3ll992l) : (8) 4Prvor & Mevlarfll993l) . 

NOTE. — The electronic edition of Astrophysical Journal Supplement Series also includes a machine-readable version of Table |13| Key to columns 



Column (1 
Column (2 
Column (3 
Column (4 
Column (5 
Column (6 
Column (7 
Column (8 
Column (9; 

Column (10)— Wilson-model A(log T v ) = log(Tj, y 7TP up ). 
Column (11) — Dynamical mass-to-light ratio calculated within power-law model fit to cluster surface-brightness profile. 
Column (12)— Power-law model A(log Tv) = log(Tj yn /TP op ). 



Cluster name. 

— Population-synthesis V-band mass-to-light ratio, from Column (6) of Table|81 
Observed stellar velocity dispersion from literature. 

Radius of effective circular aperture corresponding to velocity-dispersion measurement in Column (3). 
Source of cr p Q ^ and R ap estimates (see below). 
— Aperture-corrected line-of-sight velocity dispersion at cluster center (mean of separate results using King, Wilson, and power-law structural models). 
— Dynamical mass-to-light ratio calculated within King-model fit to cluster surface -brightness profile. 
-King-model A(log Tv) = log(Tj yn /TP op ). 

-Dynamical mass-to-light ratio calculated within Wilson-model fit to cluster surface-brightness profile. 

dyn I 
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TABLE 14 

Galactocentric Radii and k-Space Parameters 



Cluster 
(1) 



R gc [kpc] model 

(2) (3) 



(4) 



(5) 



(6) 



FORNAX 1 



MW-AM1 



LMC-HODGE 1 1 



SMC-KRON3 



123.28 



2.16 



3.71 



1.59 



K 
W 
PL 
K 
W 
PL 
K 
W 
PL 
K 
W 
PL 




NOTE. — A machine-readable version of Table ll4j is published in its entirety in the electronic edition of the Aslrophysical Journal Supplement Series; only a portion is shown here, 
for guidance regarding its form and content. Key to columns: 
Column (1) — Cluster name. 

Column (2) — Distance from center of parent galaxy, in kpc. For LMC, SMC, and Fornax -dwarf cluster s, Rg C is a projected distance (with galaxy centers given in eq. [6] of the text); 
for Milky Way globulars, R^ is the three-dimensional Galactocentric distance, taken from Harris ' 1996). 
Column (3)— Model fit. 

Column (4) — K[ = (log cr 2 + log R h )j^/2, for c p .o from Table ll2l in km s _1 and R h from Table llll in kpc. 

Column (5)— k 2 = (log cr 2 +2 log £/, - log R h )f V5, for Op.o from Tabletl2lin km s" 1 , from TabletHlin Af ^ pc~ 2 , andfy, from TablefTTlin kpc. 
Column (6) — K3 = (log er 2 - log X/, - log Rf,)/y/3, for cr p o from Table ll2j in km s" 1 , from Table ll2| in Mq pc~ 2 , and R/, from Table ll lj in kpc. 
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FIG. 12. — Fits of Kina 1 1966) models, isotropic Wilson 1 1975) models, and power-law models (eq. 0) to the surface-brightness profiles of 53 LMC clusters 
and 10 SMC clusters spanning a wide range of ages; 5 globular clusters in the Fornax dwarf spheroidal; and 84 Galactic globular clusters besides lu Centauri. 
In all panels, solid lines correspond to the Wilson-model fits; dashed lines, to King-model fi ts; a nd dotted lines, to power-law fits. All SB data are plotted with 
corrections for zeropoint changes and V-band extinction included (Columns 2 and 3 of Table llOt . R adia l scales are in arcsec along the lower horizo ntal axe s of 
all panels, and in pc along the upper horizontal axes. Numerical details of every fit are given in Table llOl and derived cluster parameters are in TableslTTIandm 
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